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We present a comprehensive study of two of the most experimentally relevant extensions of Ki- 
taev's spinless model of a ID p-wave superconductor: those involving (i) longer range hopping and 
superconductivity and (ii) inhomogeneous potentials. We commence with a pedagogical review of 
the spinless model and, as a means of characterizing topological phases exhibited by the systems 
studied here, we introduce bulk topological invariants as well as those derived from an explicit con- 
sideration of boundary modes. In time-reversal invariant systems, we find that the longer range 
hopping leads to topological phases characterized by multiple Majorana modes. In particular, we 
investigate a spin model, which respects a duality and maps to a fermionic model with multiple 
Majorana modes; we highlight the connection between these topological phases and the broken 
symmetry phases in the original spin model. In the presence of time-reversal symmetry breaking 
terms, we show that the topological phase diagram is characterized by an extended gapless regime. 
For the case of inhomogeneous potentials, we explore phase diagrams of periodic, quasiperiodic, and 
disordered systems. We present a detailed mapping between normal state localization properties 
of such systems and the topological phases of the corresponding superconducting systems. This 
powerful tool allows us to leverage the analyses of Hofstadter's butterfly and the vast literature 
on Anderson localization to the question of Majorana modes in superconducting quasiperiodic and 
disordered systems, respectively. We briefly touch upon the synergistic effects that can be expected 
in cases where long-range hopping and disorder are both present. 

PACS numbers: 03.65.Vf, 71.10.Pm 

I. INTRODUCTION 

The recent explosion of studies concerning Majorana fermions in solid state systems has brought the one-dimensional 
spinless p-wave paired superconducting wire into the limelight. As the prototype for hosting topological phases 
characterized by bound Majorana states at the ends of the wire, this superconducting system has been theoretically 
studied from a variety of angles [l|-l37j and has formed the basis for several experimental realizations [38l - l4l| . While 
these analyses have led to revisiting several theoretical aspects investigated in the past decade, such as topological 
features and symmetry classification of the system, localization properties in the presence of disorder, and physics in 
the presence of multiple channels, the prospect of experimental realization has instigated new exploration of physical 
realizations. A number of studies have focused on novel materials and geometries, as well as subjecting the system to 
controlled external potentials. Additionally, new means of detecting Majorana modes and signatures of related non- 
Abclian statistics have been proposed, as well as on applications, such as schemes for topological quantum computation 
in these systems [12, El] ■ 

Here, contributing to this vast literature, wc present a comprehensive study of the superconducting wire subject 
to various experimentally relevant modifications of Kitacv's original model. By exploring several variants of the 
^ coupling/hopping amplitudes and spatially varying electronic potentials, we build on and unify previously studied 
\ features of similar models. We investigate physics stemming from long-range hopping wherein conduction electrons 
in a lattice version of the wire can hop across several sites. In systems obeying time-reversal symmetry, we illustrate 
how such hopping may give rise to multiple end Majorana fermions. We perform an involved study of the wire in 
the presence of various potential landscapes, in particular, periodic potentials, quasiperiodic potentials and disorder. 
Finally, we briefly describe some of the richness associated with systems exhibiting both long-range hopping and 
inhomogeneous potentials. Our goal here is to present several new results of experimental relevance as well as to 
gather foundational information that is scattered in the literature, recasting some of it in simpler language and 
providing the non-expert a pedagogical, self-contained exposition leading up to the results. 

We begin with a review of the lattice version of the one-dimensional spinless p-wave superconducting fermionic 
system, namely the Kitaev chain, and focus on its Majorana mode properties. Specifically, we consider the phase 
diagram and topological features of the system by virtue of the presence of Majorana end modes in some regions of 
parameter space (topological phases) versus their absence in others (non-topological phase). Wc develop the formal- 
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ism for several alternate topological invariants that are useful in different circumstances. The relevant topological 
invariants can either be derived from bulk properties for a homogeneous system with periodic boundary conditions or 
boundary properties in long finite-sized wires which pinpoint the existence of end modes; the latter was developed in 
our previous work in Ref. |3l| within a transfer matrix formalism and is easy to apply to a range of situations. This 
connection provides an elementary illustration of the so-called bulk-boundary correspondence |44l. l45| . The topological 
invariants (TIs) are either Z2 in nature, detecting odd versus even number of Majorana end modes (and associated 
Z2 bulk physics), or Z, detecting the presence of multiple, independent Majorana end modes (and associated winding 
numbers in the bulk), depending on whether time- reversal symmetry (TRS) is obeyed (43 - l46| . It has been long known 
that the superconducting system can be mapped to a spin chain, the topological phase being associated with an Ising 
ferromagnet and the non-topological phase with a paramagnet jll.[3ll[33j: we recapitulate this mapping. 

Detailed studies involving symmetry classifications of superconductors have identified the one at hand as having a 
topological invariant that lies in Z when time-reversal symmetry is preserved (class BDI). Here, we show that long- 
range hopping explicitly brings out this Z character Depending on the strength of the hopping terms, the system 
can exhibit a slew of phases characterized by the presence of multiple, robust, independent Majorana end modes. In 
terms of bulk properties, these phases are characterized by winding numbers having different integer values. While in 
terms of lattice connectivity, the basic structure of these Majorana modes is the same as previously discussed modes 
in multichannel wires l47l l48j , our construction explicitly discusses the physics of long-range hopping and can be 
realized within a single-channel system. 

As a specific instance of a long-range hopping, we present a model having second nearest-neighbor hopping that 
also displays interesting physics in terms of spin variables. The model exhibits four distinct phases, three of which 
also have analogs in the Kitaev chain, and the fourth phase, being the novel one, supports two independent Majorana 
modes at each end of a superconducting wire. In the spin language, the model has cubic terms and an elegant 
duality map between one set of spin variables and another. The four phases correspond to long-range ordering of 
different spin components in the two sets of spin variables. The models shows a rich phase diagram highlighted by 
the multiple-Majorana phase in the fermion language and by ordering involving duality maps in the spin language. 

The symmetry classification scheme has shown that in the absence of time-reversal symmetry, the topological 
invariant associated with the system lies in Z2 (class D) j43 - l47i |49| - |5]| . Here, by explicitly describing the system 
in terms of Majorana fermion degrees of freedom, we discuss the manner in which the Z form of the invariant for 
time-reversal symmetric systems gets reduced to Z2 once symmetry breaking terms are introduced. We consider 
the effect of such terms on the Kitaev chain phase diagram, in particular, the presence of a complex phase in the 
superconducting order parameter. We show that such a phase gives rise to an unusual extended, gapless, capsule-like 
region that lies between the gapped topological and non-topological regions. 

The presence of spatially varying potentials too causes dramatic changes to the topological phase diagram of the 
superconducting wire. Expanding on our results of Ref. (soj . we develop our transfer matrix formalism and show 
that the manner in which Majorana end mode wave functions decay into the bulk of the system can be connected to 
localization properties of a normal system (vanishing superconducting gap) described by the same potential landscape 
(similar methods have been employed in 0, H^ill^)- Moreover, for a common class of potentials, we show that there 
exists a mapping between the phase boundary for fixed gap strength and its inverse. Armed with the transfer 
matrix tool, we analyze the effect of several different potential landscapes. As the simplest case, periodic potentials 
significantly affect the topological phase boundary, providing a knob to control the extent of the topological regime. 
For quasiperiodic potentials, the map to normal systems provides significant insight and reveals that the topological 
phase diagram mirrors fractal-like structures, like the Hofstadter butterfiy, that naturally emerge in normal systems 
possessing quasiperiodicity. 

Finally, for disordered potentials, the mapping to normal state properties is powerful in that it allows us to leverage 
the extensive literature on Anderson localization in normal systems to identify the topological phase diagram for 
the disordered superconductor. We consider a variety of potentials, including uniformly- and Lorentzian-distributcd 
disorder. All examples show that the topological phase continues to occupy a significant region of the phase diagram 
in the presence of disorder. One of our findings is a ubiquitous singularity in the phase boundary at the random-field 
transverse Ising critical point. While we present a fairly extensive set of phase diagrams and analyses of Majorana 
physics, our main contribution, the map to normal systems, is much more far-reaching in extent. 

Our presentation is as follows. In Sec. I, we introduce our framework in the context of a review of the Kitaev chain. 
In Sec. II, we describe the generic long-range hopping model and Sec. HI the specific instance of the four-phase 
model. In Sec. IV, we discuss the effect of time-reversal symmetry breaking. In Sec. V, in the context of spatially 
varying potentials, we further develop the transfer matrix technique introduced in Sec. I and apply the methods to 
the case of periodic and quasiperiodic potentials. We analyze disordered potentials in Sec. VI and conclude in Sec. 
VII. 
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II. REVIEW: TOPOLOGICAL ASPECTS OF A p-WAVE SUPERCONDUCTING WIRE 

We review the salient features of a single-channel p-wave paired superconducting wire within the context of the 
broadly- used ID tight-binding system of spinless electrons pioneered by Kitaev [l| (referred to as Kitaev chain). We 
approach this simple and well-studied system from various angles as a preparation for the new material in subsequent 
sections. We discuss the phase diagram of the Kitaev chain in terms of its topological properties characterized by the 
presence or absence of zero energy Majorana modes at the ends of a long and open chain. To establish these topological 
properties, we consider two symmetry classes, one respecting time-reversal invariance and the other breaking it. To 
analyze these classes, we present various topological invariants, some exploiting bulk features of the Hamiltonian in 
momentum space and others explicitly counting the number of zero energy modes at the ends of an open chain. 
Finally, we review the Jordan- Wigner transformation [s^ which maps the fermionic system to a spin- 1/2 chain. 



A. Kitaev Chain: Model, Dispersion and Phases 



In this lattice description of the single-channel p-wave superconductor, electrons experience a nearest-neighbor 
hopping amplitude w, a superconducting gap function for pairing between neighboring sites. A, and an on-site 
chemical potential /i; in this section, we will assume that all these parameters are real. For a finite and open wire 
with A/" sites, the Hamiltonian takes the form 



H = 



E 



AT 

{Pnfn+1 + fl+Jn) + A(/„/„+i + fl^j}, ) ] - E (/"/« " 1/2) , 



(1) 



where the operators /„ satisfy the usual anticommutation relations {/,«, fn} = and {/,„, /t} = S„in- We can assume 
that w > 0; if w < 0, we can change its sign by the unitary transformation /„ — > (— 1)"/„. (Throughout this paper 
we will set both h and the lattice spacing equal to unity). It is important to note that Majorana end modes can 
only appear if the superconducting order parameter A 7^ 0. This is because Majorana modes do not have a definite 



fermion number, while the Hamiltonian commutes with the total fermion number, /n/m if A = 

Towards exploring the Majorana mode structure of the wire, we can decompose the electron operator in terms of 
real (Majorana) operators a and b as 



fn 



(2) 



The operators a„ and 6„ are Hermitian and satisfy {a,„, a„} ~ {&,„, b„} ~ 2(5„i„ and {a„i, 6„} = 0. Then Eq. ([T]) can 
be re-written as 



-ff = ^ [(it; - A)a„5„+i + (it; + A)a„+i6,; 



n=l 



AT 
n=l 



(3) 



To study the bulk features of the system, we consider a long wire having periodic boundary conditions (so that the 
first summation in Eq. ([l]) goes from n = 1 to Af). Then the momentum fc is a good quantum number and it goes 
from — TT to TT in steps of 2-k jM . Defining the Fourier transform /fc ~ Yl^=i fn e**^", Eq. ([T]) can be re-written 
in momentum space as 



H = 



E ( fl ) ( /t 



0<fe<7r 



fk 



'(2wcosfc + /i) -I- 2Asinfcs^, 



(4) 



where the s° are Pauli matrices denoting pseudo-spin degrees of freedom formed by the fermion particle-hole subspace. 
The dispersion relation follows from this and is given by 



Ek = ± 



\J (2wcosfc + ju)^ 



4A2 sin^ k. 



(5) 



The energy vanishes at certain values of fc; these are given by lines in the two-dimensional space of the parameters 
and A/w. These gapless lines correspond to phase transition lines which separate different phases. The phase 
diagram consists of three lines demarcating phases I, II and HI, as shown in the top left diagram in Fig. 21 On the 
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vertical red lines lying along /i/w = ±2, the energy vanishes at fc = tt and zero respectively, while on the horizontal 
blue line extending from fi/w = —2 to 2 at A/w = 0, the energy vanishes at fc = cos~^{—iJ,/{2'w)). 

Phase Diagram- To obtain insight into the nature of the phases, we can consider some extreme limits. As shown 
in Fig. [T] (a), for A = w = 0, ^ ^ 0, which lies in phase III of the phase diagram, we see in terms of the Majorana 
mode Hamiltonian of Eq. ([3]) that each Majorana mode a„ on a given site is bound to its partner 6„ with strength 
/i, leaving no unbound modes. For A = w 7^ 0, /.j = 0, lying in phase I, the only existing bonds connect a„ to its 
neighbor bn+i (Fig. [1] (b)), leaving a free a/6-Majorana mode at the right/left end of a finite sized system (Fig. [1] 
(c)). For A ~ ^ 0, /Li = 0, lying in phase II, the roles of a and b modes become interchanged. As shown using 
topological arguments in the next section, the presence/absence of these end modes is robust in that deviations from 
these extreme limits in parameter space does not change these features unless a phase boundary associated with a 
vanishing gap is crossed. 



a 
b 
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(d) (e) 

FIG. 1: Pictures of a fermionic chain showing the couplings between Majorana operators a„ and 6„ at different sites when only 
one kind of interaction, ambn, is present in the Hamiltonians in Eqs. ((3]) and (|16[) . Panels (a), (b) and (d) respectively show 
couplings of the form a„b„, a„6„+i (Eq. ((Sj) and a longer range interaction a„b„+2 (Eq. (|16p ) in the middle of a long chain. 
Panels (c) and (e) show couplings of the form a„6„-|-i and a„6n+2 for a chain which has been cut into two, thus giving two open 
chains. Panel (c) has a Majorana mode 03 at the right end and a mode 64 at the left end of an open chain. Panel (e) has two 
Majorana modes as and 04 at the right end of an open chain and two modes 65 and be at the left end of an open chain. (Panel 
(a) would not have any Majorana end modes for a chain cut into two). 



B. Topological Invariants (TIs), Transfer Matrix Approach and Majorana Modes 

The topological properties of the superconducting wire described above can be captured in a variety of ways, all 
involving global features of the system characterized by TIs. There are several TIs developed in the literature; some 
are a type of generalized winding number (for instance, the celebrated TKKN invariant of the integer quantum Hall 
effect [Hj] ) or are based on a Pfaffian, as in [ij . Building on previous work [13, HH , we introduce complementary 
'boundary invariants', i.e. TIs which are derived from an explicit counting of the number of edge Majorana states. 
It should be emphasized that although the bulk and boundary invariants are derived in different ways, they all share 
two crucial characteristics: (1) they are restricted to integer values and (2) small deformations of the Hamiltonian 
which do not close a bulk gap cannot change their values. Furthermore, they enumerate the number of edge Majorana 
modes. The different forms will be used in subsequent sections based on convenience and the aspects studied. 

Symmetry Classes- In this paper wc will encounter both 1 and Z2 TIs (which we explain below). These 
types indicate the level of topological protection enjoyed by a topological insulator or superconductor. An extensive 
classification of topological insulators and superconductors has been developed, and the type of TI can be determined 
from the class which the Hamiltonian falls into [43 - l46| . The class depends on a number of factors, such as whether 
the particles are bosons or fermions and whether the Hamiltonian has spin rotation symmetry (if the particles have 
a non-zero spin), particle-hole symmetry or TRS. The different symmetry classes have been enumerated and there is 
extensive discussion in the literature [47l. l49l - [5ll [54| . Relevant to this work are the classes: 

Class D - This class is ap prop riate for describing spinless electrons exhibiting p-wave pairing but for which time- 
reversal symmetry is broken |47l. I49l - l5l1 | . This class exhibits Z2 topological protection (in ID) indicating that the 
parity of the number of Majorana modes is protected p]. |43 - I46j . 
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Class BDI - This class also describes spinless electrons exhibiting p-wave pairing, with the additional restriction 
that time-reversal symmetry is obeyed. This class is associated with a Z TI in ID, indicating that the number of 
Majorana modes itself is topologically protected (47l |. 

As is seen here, a general feature of the 'periodic table of topological insulators' is that the level of topological 
protection is generally greater as the symmetry constraints on the Hamiltonian become more restrictive [U113- 



1. Bulk Invariants 



Depending on the time-reversal properties of the Hamiltonian h^, two different invariants can be defined. (Although 
Eq. Q defines hk only for < < tt, it will now be convenient to use the same expression for hk for the entire range 
< k < 27r). In Eq. (j?]), we see that hk is of the form hk = a^kS^ + o,2kS^ , which maps to the vector Vk ~ a^kZ + a2ky 
in the z — y plane. The form is general to this class of systems which have TRS, i.e., when h*_^, = hk for all k. 
Our first bulk invariant W ^ as introduced in Refs. fssl [55| . is then the winding number associated with the angle 
(j)k ~ tan^^(a2fe/a3fe) made by the vector Vk with respect to the z axis upon traversing the Brillouin zone, i.e., 

w ^ r ^ ^. (6) 

Jo dk ^ ' 

This object can take any integer value and is a TI, namely, it does not change under small changes in hk unless hk 
happens to pass through zero for some value of fc; in the latter case, the winding number becomes ill-defined and the 
energy Ek = which means that the bulk gap closes at that value of A:. As will be explicitly discussed for the case of 
long-range hopping, we thus have a Z-valued invariant, as expected from the general classification discussed above. It 
is straightforward to determine the relationship between W and the number of edge modes Na and Nf, of a and b type 
(respectively) of Majorana modes at the left-hand end of an open chain, and confirm that W = Nt, — Na- Note that 
taking A — > —A reverses the winding of Vk and thus W, which is consistent with this transformation interchanging 
a and b Majoranas (see Sec. Ill A| ). 

In the absence of TRS, i.e., if /i^j, ^ hk for some fc, as seen in later sections, hk generally has four components, 
hk = aokl + aikS^ + a2kS^ + o-skS^, and it is not possible to define a winding number for the corresponding vector 
in four dimensions. On the other hand, for momentum values k = and tt, hk has only one component, namely, 
ho = h{0)s^ and h^^ ~ h{TT)s^. As seen in Eq. (|4]), this is naturally true also in the case with TRS. Under the stringent 
assumption that the system is fully gapped (i.e., gapped for all values of fc), one can define a Z2-valued TI: 

i^buik = sgn (/i(0)/i(7r)) . (7) 

This invariant can take the values —1 or 1 and is topological in that its value cannot change unless either h{0) or h{Tr) 
crosses zero in which case the energy Eq or i?^ vanishes. We will see that Vbuik is equal to the parity of the number 
of Majorana modes at the end of the system. 

The fact that systems which break TRS have a symmetry only at the level Z2 may also be seen another way. As 
seen above, sending A — >■ —A takes W — >■ —W. This process can be carried out for class D without closing a gap 
since we can take A = Aoe*"^ and take cp from — tt. However, for class BDI, A is constrained to be real and the 
only way to take A — A is to pass through zero, which closes a gap (47j . 



2. Boundary Invariants: Transfer Matrix Approach 

We now outline the transfer matrix approach detailed in our previous work, Ref. [sH ] , for identifying the Majorana 
modes structure at the end of a wire. It will become apparent that the transfer matrix explicitly gives us Na and 
Nb (as defined in the previous section, Na^b is the number of a- and &-type Majorana modes at the left-hand side of 
the system, respectively) . We also note that finding Na and Nb from the transfer matrix will only involve finding the 
eigenvalues of a matrix; this is numerically easier than calculating than W which involves doing an integral over the 
momentum k. 

The transfer matrix can be obtained from the Heisenberg equations of motion for the Majorana operators in Eq. ([3|) : 

[iJ, a„] = -i(w + A)6„+i - i(tii - A)6„_i - i/^&„, 

[H,bn] = i(w' + A)a„_i + i{w - A)an+i + iiian- (8) 



. dbn 
^~dt 
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Assuming that these operators depend on time as a„ = a„e~*^* and 6„ = /?jie~'^*, we find the values of the energy 
E for which the above equations have solutions. The solutions in the bulk have the same dispersion as the one given 
in Eq. ([5]); in particular, the energies difi^er from zero by a finite gap except on the phase boundaries. In addition 
to these bulk modes, in the topological phases there are end modes which lie at zero energy for long chains and are 
therefore separated by a gap from the bulk modes. (For short chains, modes at the two ends can hybridize and lift 
their degeneracy away from zero.) For E ^ 0, Eqs. ([8]) take the form 



{w + A)a„+i + {w - A)a„_i + = 0, 

(w + A)^„_i + (w~A)/3„+i + ^il3n = 0. 

These equations can be represented in the transfer matrix form 



an 



where A„ = 



At A— til 
" A+w A+to 



1 







(9) 



(10) 



Since the An may be taken as functions oi n/w and A/w, we set w = 1. A similar expression holds for the transfer 



matrix i?„ for the /3„ 
and Bn vary with n). 



(In a later section, we will consider models where the chemical potential /i„ and therefore A„ 




FIG. 2: Panels (a-c): Possible locations of the eigenvalues of a 2 x 2 transfer matrix for Majorana mode a as discussed in Sec. 
II B. Only one eigenvalue may lie inside the unit circle as in (a), or both eigenvalues may lie inside the unit circle as in (b), 
or both may lie outside the unit circle as in (c). Panel (d-g); Possible locations of the eigenvalues of a 3 x 3 transfer matrix 
for Majorana mode h as analyzed in Sec. V A. As discussed in the text, the number and type of Majorana modes at each end 
of a long chain is governed by the number of eigenvalues inside the unit circle; cases (a) and (e) give non-topological phases 
characterized by the absence of end modes, while all the other cases give topological phases with one or more end modes at 
each end. 



The existence of end Majorana modes requires the a„ (or /?„) to be normalizable, i.e., ^„ |a„p (or ^„ |/3nP) 
should be finite. The number of eigenvalues of the full transfer matrix A = Il^i with magnitude less than I is 
denoted by ny. The number nf corresponds to the number of roots of the characteristic polynomial for the transfer 
matrix, f{z) ~ det (/ — Az), that lie within the unit circle. Hence, 

""^"i/l^/^^'^^^/^^^^' ^^^^ 

as was noted in the context of plane wave zero modes in Ref. [s^. One can show that the eigenvalues of the full 
transfer matrix B = Jl^i are inverses of the eigenvalues of A, and therefore the number of eigenvalues of B with 
magnitude smaller than 1 is 2 — nf. For nj = and 2, q;„ and /3„ are normalizable and therefore the system is 
topological (with a a mode at one end and a /? mode at the other end of a long chain), whereas for nj ~ 1, a„ and 
/3„ are not normalizable and the system is non-topological. 
With nf at hand, we can define a Z2 invariant 



^end = - (-1)"^ =-Sgn(/(l)/(-l)), 



(12) 
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for which Vend = — 1 and Vend = 1 reflect an odd versus even number of normalizable end Majorana modes, respectively. 
For the Kitaev chain of Eq. ([T]), since there can only be one or zero modes, these values correspond to topological 
and non-topological phases, respectively. Since the topology of the system depends only on the magnitude of A, we 
take A to be positive; hence | det „4| < 1. Then the two eigenvalues of A obey IA1A2I < 1. Therefore, for |Ai| < IA2I, 
we have |Ai| < 1 and rif is completely determined by the larger eigenvalue A2. Thus, we have that v = sgn(ln|A2|) 
for A > 0. We claim that the odd-even Majorana structure denoted by the Z2 invariant Vend of Eq. p2p matches the 
form of the bulk invariant Vbuik of Eq. ([7]). 

Phase Diagram - Wc can now revisit the Kitaev chain phase diagram in light of these invariants. Considering 
some special cases, for /i = and A > 0, we find that both eigenvalues of A have magnitude smaller than 1, so that 
Uf = 2. Further, we find that there is a zero energy Majorana mode of type a at the left end and of type b at the 
right end of a long chain. This describes phase I following the discussion after Eq. ([5]). If jU = and A < 0, both 
eigenvalues of A have magnitude larger than 1, so that Uf = 0. We find that there is a zero energy Majorana mode of 
type b at the left end and of type a at the right end of a long chain. This describes phase II. Finally, let us consider 
A = and > 2t. This lies in phase III. We then find that n/ = and that there are no Majorana modes at either 
end of a long chain. This is summarized in Table I. 



Phase 


1/ 


W 


nj 




Nb 


I 


-1 


-1 


2 


1 





II 


-1 


1 








1 


III 


1 





1 









TABLE I: Summary of the topological phases (I and II) and non-topological phase (III) exhibited by the Hamiltonian of Eq. HI) 
and the different topological invariants. Here u refers to both Vbuik and Uend, consistent with their equivalence. 



Note that on the phase transition lines which separate the various phases, the bulk gap closes, i.e., Ek = at 
certain real values of the momentum k. Eqs. (|9]) and (jlOp then imply that An must have an eigenvalue of the form 
e^'' which lies on the unit circle. Thus the eigenvalues of the transfer matrix cross the unit circle as we go across a 
phase transition line. 

The TIs that we have discussed here are related to other invariants which exist in the literature. We refer the reader 

to Refs. [Uliilii. 



C. Mapping to a Spin-1/2 Chain 



The Hamiltonian in Eq. (|3]) can be mapped to that of a spin-1/2 chain; here we briefly review the mapping and 
associated physics [s^]. We define the Jordan- Wigner (JW) transformation between a spin-1/2 and a spinless fermion 
at each site n so that the states with — ±1 correspond to the fermion number = 1 and respectively. The 
JW transformation takes the form 



n 



(-1)" 



(13) 



Eq. ^ can then be re-written as 



Af-l ^ N 



n=l 



(14) 



with 



= (w - A) /2, Jy = {w + A) /2, h = pL. 



(15) 



The Hamiltonian describes a spin chain having nearest-neighbor xx and yy couplings and a magnetic field pointing 
in the z direction. 

In the discussion below Eq. ([5]), we had stated that there arc three phases, I, II and HI, depending on the values 
of the parameters ^/w and A/iu. While in the fermionic language the ordering is topological in nature and has no 
local order, we can see that in the spin language, the corresponding phases are described by ordering of local spin 
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variables. Assuming that w > 0, let us first consider ^ = 0. Then the point A = w lies in phase I and corresponds 
to Jj; = and Jy = w > 0. This describes a ferromagnetic Ising model with yy couplings which has long-range order 
in a-y . Similarly, A = lies in phase II and describes a ferromagnetic Ising model with xx couplings which has 
long-range order in a^. Finally, A = and \fi\ ^ 2w describes a model with = Jy and a magnetic field in the 
z direction whose magnitude is much larger than J^. This describes a system with no long-range order in either cr^ 
or . These special points represent the entire phase diagram in that phase III is that of a disordered paramagnet, 
and phases I and II have Ising ferromagnetic order (along y and x, respectively) and arc separated by a gaplcss line 
(A = 0, Jx — Jy) describing an xy ordered spin chain in a transverse field. 

While the spin chain described above has a venerable history in and of itself, the mapping sets the stage for more 
complicated spin systems and mappings in the context of topological order, such as with the Kitaev honeycomb model 
and the Kitaev ladder [sil, [32|. In Sec. IV, we will discuss an interesting generalization of the fermionic model, 
which has multiple Majorana modes and exhibits a rich phase diagram, that can be elegantly described in terms of 
spin variables. 



III. LONG-RANGE HOPPING AND MULTIPLE MAJORANA END MODES 



Here we argue that the presence of longer range hopping that extends beyond the nearest neighbor has the dramatic 
consequence that in a spinless superconducting wire that preserves TRS, multiple topologically protected Majorana 
modes can form at the ends of the wire. These modes are stable in that regions in parameter space corresponding to 
different number of modes are protected by bulk gaps in the energy spectrum and correspond to topologically distinct 
phases. In the previous analyses of symmetry classes mentioned above, it has been argued in various ways that the 
TI lies in Z for the class currently being considered, class BDI. The discussion here in terms of long-range hopping in 
a single-channel model and of the TIs defined above gives a simple, direct and comprehensive picture for the Z form 
by way of multiple Majorana end modes. 

Below we introduce the general long-range hopping model and analyze its topological properties based on the 
presence of multiple Majorana end modes. In the next section, we explore the features of the aforementioned model 
having four topologically distinct phases and a dual Ising representation in the language of spins. 



A. Long-range Hopping Model 



As a very general case of long-range hopping, we consider a modification to the Kitaev chain of the previous section 
that takes into account an infinite set of couplings for long-range hopping. Focusing purely on the hopping and 
pairing, this generalized version of the Majorana Hamiltonian of Eq. ([3|) for an infinite wire takes the form 

00 00 

-ff = - i ^ ^ a„6„+r, (16) 



r— — 00 n— — 00 



where the Jr are real parameters. Note that we have chosen to discuss long-range hoppings in a single chain. However, 
our model is equivalent, for appropriate choices of the J^, to a multi-chain system with both interchain and intrachain 
couplings. 

The Hamiltonian in Eq. (jl6p is invariant under TRS which involves complex conjugating all numbers, taking the 
time t — >■ —t, changing a„ — > — a„ and keeping &„ unchanged. The transformation of a„ is justified by Eq. (|13p where 
we see that complex conjugating reverses the sign of cr|J and therefore of a„. Note that the square of the time- reversal 
transformation is equal to -1-1 since it leaves a„ and 6„ unchanged. 

In terms of Dirac fermions = (a„ -|- i6„)/2 and = (a„ — i6„)/2, the generalized momentum space version of 
Eq. (|4]) momentum space stemming from Eq. p6|) takes the form 

0<k<7r 
00 

hk = -2 E [JrCOs{kr) s"" + JrSm{kr) s^], (17) 

r— — oo 

where the are once again Pauh matrices. 
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Note that the Hamiltonian is time-reversal invariant: 
Eq. (dZD, 



'-A; 



hk- The energy- momentum dispersion follows from 



Eh 



±2 



Jr cos{kr) 



(18) 



The phase diagram can, in principle, be found from this expression by demanding that should vanish for some 
value of k lying in the range [0, tt] or by employing the TIs defined in the previous section. While this is a difficult 
problem, we will see below that a great deal of insight can be gained by considering situations in which one of the Jr 
is much larger than all the others. 



B. Multiple Majorana Modes and Topological Phases 



To address what phases are exhibited by the model for an infinitely long chain and how many zero energy Majorana 
modes there are for the ends of a long chain in the different phases, consider a Hamiltonian where one of the Jr, say 
Jq for some positive integer q, is non-zero and positive while all the others are zero [47| . (From Eq. (jl6[) we see that 
the ground state is then given by a state in which ianbn+q — 1 for each n). We show this situation in Fig. [1] (d) for 
the case of g = 2. In Fig. [1] (e), we have shown a dotted line which cuts the chain into two; it is clear that the open 
chain on the right side has q Majorana modes of type h at its left end, while the open chain on the left side has q 
Majorana modes of type a at its right end. Since we expect the number of Majorana modes and their types to be 
TIs, this phase will survive for small changes in the values of all the other J^. We therefore conclude that the model 
in Eq. (|16|) has an infinite number of phases which can be labeled by an integer q which describes the number and 
types of Majorana modes at the ends of a long chain. 

Next, we look at what the Z- valued bulk invariant, W , defined in Eq. ([6]) gives for the Hamiltonian in Eq. (|17l) . 
As before, we can think of hk as defining a two-dimensional vector. If only one of the couplings Jq is non-zero and 
positive, the vector is given by 

Vk ^ -2 [Jq cos(fc(?) z + Jq sm{kq) y]. (19) 

As k goes from to 27r, this generates a closed curve which encircles the origin of the z — y plane q times in the 
clockwise direction. Defining (j)k = tan~^(T^/T4) and using Eq. we find that the winding number is equal to q for 
the configuration given in Eq. (|19p . i.e., W = q. Thus the winding number is also equal to the number of Majorana 
modes of type b (a) at the left (right) end of a long chain. We can now consider what happens if all the Jr are allowed 
to be non-zero; the vector is then given by 

Vk = -2^ [Jr cos(fcr) z + Jr sm{kr) y]. (20) 

r 

The winding number is a TI and therefore does not change under small changes in all the Jr- The system thus 
continues to remain in the phase q as long as the closed curve does not pass through the origin for any value of k. 
But if the curve passes through the origin for some value of fc, the energy vanishes at that value of k and the system 
lies at a quantum critical point separating two phases having different topological values q. 

Transfer matrix approach with constraint equations - We can also understand the existence of Majorana 
end modes using the transfer matrix approach. Consider a semi- infinite chain, with sites going from n = 1 to oo, with 
a Hamiltonian of the form 

oo g oo 

^ = E E - ^ Mn(/,t/n-l/2), (21) 

n— 1 ?n — 1 n— 1 

where Wm-, A™ and fin are all real, and g > 1 is an integer. Let us now study the zero energy equations of motion to 
see if there are Majorana modes localized near n = 1, i.e., the left end of the chain. To be specific, let us focus on the 
a„ modes. Each site n > q + I gives rise to an equation linking a„ to all sites from a„-q to Un+q] these equations can 
be described by a d x d transfer matrix A, where d = 2q. However, for n = 1, 2, • • • ,q, all the sites up to a„_q are 
not present in the system; hence the corresponding equations are not of the transfer matrix form, but instead provide 
q constraints on the first 2q values of a„. Let us now suppose that the parameters {wm-, ^m, fJ-m) are such that the 
transfer matrix A has Uf eigenvalues with magnitude smaller than 1, and the other d—Uf eigenvalues have magnitude 
larger than 1. The eigenvectors corresponding to the first Uf eigenvalues are normalizable. However, the presence 
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of q constraints means that there are only Uf — q independent and normalizable end modes. Thus, the number of a 
Major anas on the left-hand side of the system is 

Na = nf-q. (22) 

If Uf < q, there are no Majorana end modes of type a. A similar analysis carried out for the 6„ modes indeed shows 
that there are an equal number of b Majoranas localized to the opposite end of the system. Finally, we also note that 
this argument implies that a and b type Majoranas can never occur on the same side of the system. Such a state is 
incompatible with TRS. As illustrations, models of the form given in Eq. (|2ip with next-nearest neighbor hopping 
and multiple Majorana modes are presented in Fig. [T3]as well as in Sec. IV. 

While the above analysis gives the exact number of end Majorana modes of types a and 6, it can also be useful 
to find an expression for the Z2-valued invariant analogous to Eq. ((T^ . We can derive this as follows. Since the 
transfer matrix A is real, its eigenvalues must either be real or must come in complex conjugate pairs. If we define 
/(A) ~ det{A— A/), we have the relation (— 1)"^" ~ sgn (/(l)/(— 1)). (This relation holds because if an eigenvalue of 
A, denoted by A^, is a real number not equal to ±1, then Af — 1 < (> 0) depending on whether A^ is smaller than 
(larger than) 1 in magnitude. If Xi is complex, then (Af — 1)(A*^ ^ 1) > 0). If we define 

z., = (-l)^sgn(/(l)/(-l)), (23) 

we see that Vq = (—1) corresponds to having an even (odd) number of Majorana modes of type a at the left end 
of a chain. In particular, z^g = — 1 means that there is at least one Majorana end mode and therefore the system is in 
a topological phase. Thus Eq. (|23p is the appropriate generalization of i^end (Eq. (|12p ) to arbitrary values of q and Uf 
Finally, a few comments are in order on the general long-range hopping model. We have seen that when only 
couplings between a and b modes are present, the system has Z topological symmetry, allowing for an arbitrary 
number of end Majorana modes. When TRS is broken, however, intra-couplings between the a's or the 6's themselves 
become manifest and, as seen is following sections, this couples modes near each end, allowing for only a Z2 symmetry. 
On another note, as for mappings between the fermions and spin outlined in the previous section, Eq. (|2ip yields spin 
Hamiltonians typically containing multi-spin terms involving arbitrarily long strings of . While these general cases 
are too complex to provide further insight, we present a model having a tractable mapping to a spin model in the 
next section. 



IV. LONG-RANGE HOPPING: ISING DUALITY AND A FOUR-PHASE MODEL 



As an illustration of the long-range Majorana hopping model discussed in the previous section, we now make a 
detailed study of a model with four parameters and four associated phases. We present this model not only as an 
instance of supporting multiple end Majorana modes but also as a novel spin system that supports an Ising duality 
and interesting interpretations of phases in terms of spin ordering. The spin Hamiltonian that we introduce not only 
has usual linear and quadratic terms in the spin language, but also an unusual cubic term, and is of the form 



H 



(24) 



Taking the site label n to run over all integers in Eq. (fM)) . we can define a dual spin- 1/2 chain whose site labels run 
over n + 1/2, with the mappings 



'n-l/2^n+l/2 Si+3/2 



Eq. ([M]) then becomes 



H 



'ri-l/2~n+l/2 — 

V y 



[^'''n- 1/2 ''n -I- 1/2 + + /2'''n-|-l/2 Jy'''n-1 /2'''n+l /2'''n+Z/2\ 



(25) 



(26) 



which interchanges Jj. o /i and Jy ^ 1/ with respect to Eq. ([M)) : this is the duality property. 

We note here that a model with the Hamiltonian given in Eq. but with only three parameters, ^, Jx and 
Jy, was elegantly analyzed in Ref. [33|. Since that study did not include the parameter i/, it did not enjoy complete 
duality; further, the transfer matrices considered there were two-dimensional (as in our earlier sections) rather than 
three-dimensional as we will discuss below. 
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As a map to a fcrmion model, we invoke the Jordan- Wigner transformations of Eq. (jl3|) . 
a„ = (n"=-oo'''f) '''n a-nd fo„ = (lYi^-oo'^i^ (—1)" ct^^, to obtain the Hamiltonian 

H = - « ^ [i^a„6„+2 + Ja;a„6n+1 + M^n^n + JyanK-l] (27) 

In principle, we can perform a Jordan- Wigner transformation from the r variables as well, and the resultant Majorana 
Hamiltonian also makes the duality manifest. 



A. Majorana Modes, Spin Ordering and Phases 

As in the previous section, we identify zero energy modes employing Heiscnberg equations of motion for the Majorana 
operators: 

iybn+2 + JxK+l + M^n + Jybn-1 = 0, 
Jyan+1 + M^n + Jxan-1 + VCLn-2 = 0, (28) 

with similar equations for dn+1/2 f^nd Cn+1/2 respectively. Taking 6„ = A", we obtain the cubic equation 

+ J^\^ + ^(A + Jy = 0. (29) 

This equation is essentially the characteristic polynomial equation det{B — XI) = for the 3x3 transfer matrix B 
which relates (fo„+2, ^n+i, ^n) to (6„+i, 6„, 6„_i). 

This polynomial equation has three roots, at least one of which must be real, and the complex roots must come in 
complex conjugate pairs. (The corresponding equation for a„ = k" is given by JyK^ + fiK^ + JxK + ly = Q. The roots 
of this are clearly the inverses of the roots of Eq. (P^ . and it suffices to know the roots of one equation to deduce that 
of the other). Note that if all the four parameters in Eq. pQ]). v, Jx, Jy, are scaled by the same factor, the roots 
remain the same. Further, just as the three roots are completely determined by the values of the four parameters, 
the four parameters are also completely determined, up to an overall scale factor, by the value of the three roots: if 
the three roots are Ai, A2, A3, we must have 

(A - Ai) (A ~ A2) (A - A3) = A^ + -A2 + ^X + :^. (30) 

Assuming that none of the roots A^ lie on the unit circle, there are four possible cases as shown in Fig. [5] (d-g). 
By considering certain extreme limits of the couplings, based on which spin ordering dominates and on the transfer 
matrix structure for the a and h modes, we can characterize the four regions as follows. 

Phase A- All three roots lie inside the unit circle. This is the case if \v\ ^ | J^|, \Jy\. Eq. ([25]) then implies 
that the Z2 symmetry t^_|_]^/2 ~^ ~''"n+i/2 spontaneously broken and T^j^ij^ develops long-range order. (If — >■ cx), 
we have -|-1 or —1) at all sites, while xiv ^ —00, we have an a ferromagnetic state in which t^^^^2 same value 

(either antiferromagnetic state in which t^^^^2 takes the values -1-1 and —1 on alternate sites. In either case, t^^_^_i/2 
has long-range order). As for the end modes, we can use the arguments based on the transfer matrix approach with 
constraint equations in Sec. IV B to show that this phase has two Majorana modes of type b at the left end and two 
of type a at the right end. 

Phase B- Two of the roots lie inside the unit circle while one lies outside. This is the case if \ Jx\ ^ |z^|, | 
Eq. ([M)) then implies that the Z2 symmetry cr^ — >■ —a^ is spontaneously broken and af-^ develops long-range order. 
Similar arguments as above show that this phase has one Majorana mode of type b at the left end and one of type a 
at the right end. 

Phase C- One of the roots lies inside the unit circle while two lie outside. This is the case if ^ \Jx\, W\, \Jy\- 
Eq. (j26|) then implies that the Z2 symmetry tI^^^/2 ^ spontaneously broken and t^^-^^^^ develops long-range 

order. This phase has no end Majorana modes and is therefore non-topological. 

Phase D- All three roots lie outside the unit circle. This is the case if | Jy | ^ \Jx\, W\- Eq- ([Ml) then implies 
that the Z2 symmetry tT-( —a^ is spontaneously broken and develops long-range order. Here, the end mode 
structure is switched compared with case (B) in that there is one Majorana mode of type a at the left end and one 
of type b at the right end. 

We remark that this model, which has a duality property in the spin values, is not symmetric under an exchange 
of a and b Majorana modes; the Hamiltonian in Eq. (|27l) lacks a term of the form a„6„_2. As a result, there is no 
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phase analogous to phase A in which there are two Majorana modes of type a at the left end and two of type b at the 
right end of a chain. 

We now discuss the stability of the phases under changes in the different parameters. Note that all points (Ai, A2, A3) 
in region A can be smoothly taken to each other without any of the roots crossing the unit circle, subject to the 
restriction that at least one of them is real and the complex roots come in conjugate pairs. (Note that a root can be 
made to pass through zero from the negative real side to the positive real side by taking Jy through zero. Similarly, 
a root can be made to pass through 00 from the very large positive real side to the very large negative real side by 
taking v through zero. In either case, the root does not pass through the unit circle). Since the three roots and the 
four parameters (up to a common scale) are related to each other by Eq. (pO|) . we see that it is not necessary that a 
point in region A must correspond to v being much larger than the other three parameters. However, any point in 
region A can be smoothly taken, without crossing a phase boundary, to a region in which v is much larger than the 
other three parameters. 

The unit circle |A| = 1 corresponds to real values of fc in A = e**^, i.e., the energy vanishes at some real momentum 
k lying in the range [0,27r]. So the condition that none of the roots lie on the unit circle means that the energy is 
gapped away from zero. The fact that all points in region A are smoothly connected to each other without crossing 
the unit circle means that they are all in the same phase, i.e., t^^^^^ has long-range order. On the other hand, to go 
from one phase to another (say, from A to B), at least one of the roots must go through the unit circle, so that the 
energy must at some point touch zero for some real momentum. 

To sum up, we have shown that the model at hand has four distinct phases characterized by spin ordering of a^^^ 
or T^/^ operators. In terms of Majorana end mode structure, phase (A) is the most interesting in that it is the only 
one distinct from those found in the simple Kitaev chain and it hosts two independent Majorana modes at each end. 



B. Phase Diagram 



Having identified the phases, we now study the phase diagram for different values in parameter space. Since the 
phase does not change if all the four parameters are scaled by the same number, we only have a three-dimensional 
parameter space to consider. For convenience, we recast the four couplings in terms of new parameters: 

Jx = J cos 01, Jy = J sin 01, 

^ = M COS02, ly ~ M sin02, (31) 

enabling us to study the phases as functions of 01, 02 a-nd J/M. We present the resultant phase diagrams in Fig. [3] 
where the nine panels correspond to different values of J/M . In each figure, the x and y axis correspond respectively 
to 01 and 02 lying in the range to 27r. Note that every figure is invariant under 0i — >■ 0i -I- tt or 02 — > 02 + tt or 
both. This is because these transformations either do not change any of the roots or change the signs of all the roots, 
as we can see from Eq. ([29|) and both situations leave the phase unchanged. 

We see from Fig. [3] that only phases B (green) and D (yellow) survive for J/M — > 00, and only phases A (blue) 
and C (red) survive for J/M 0. There is a critical value of J/M = (1 -I- •\/5)/2 ~ 1.618 above which phase A is 
completely absent. This is given by the fact that if the three roots Ai, A2 and A3 all lie on or within the unit circle, 
then one can show that 

J_ ^ / (Ai + A2 + A3)^ + (AiA2A3)^ 

M Y l + (AiA2 + A2A3 + A3Ai)2 ^ ^ 

has a maximum value of (1 + V5)/2 corresponding to (Ai,A2,A3) = (1,1,— (3 — V5)/2) or (—1,-1,(3 — V5)/2). 
Similarly, there is a critical value of J/M = (-s/S — l)/2 ~ 0.618 below which phase D is completely absent. Hence 
the golden ratio determines the critical value of J/M for hosting more than three phases. 

The most complex figure corresponds to J/M — 1. We see that the phase boundary lines correspond to straight 
lines crossing each other at 90°. This can be understood analytically as follows. On a phase boundary, we have 
A = e**^, with real k, in Eq. (|29p . li k ^ or it, we can take the real and imaginary parts of the equation to obtain 

{v + Jy) cos(3/s/2) + {ji + Jx) cos(fc/2) = 0, 

{v-Jy) sin(3fc/2) - (/_i-J^) sin(fc/2) = 0. (33) 

Using the identities cos(3fc/2)/ cos(fc/2) = 2cosfc — 1 and sin(3fc/2)/ sin(fc/2) ~ 2cosfc -f- 1, we obtain cosfc ~ [v + 
Jy ^ Jx ^ ^j) / {^{i' + Jy)) = — Jx + Jy — i^) / {2{v — Jy)). Eliminating cos A: from the above equations, we find 
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J/M = 5 J/M = 1.6 J/M = 1.3 




A B C D 



FIG. 3: Phase diagrams of Eq. {241 for J/M = 5, 1.6, 1.3, 1.1, 1, 1/1.1, 1/1.3, 1/1.6 and 1/5. The different colors 
represent phases characterized by different number of Majorana modes at each end of a finite system in the fermion language 
and long-range ordering of different operators in the spin language. These phases correspond to A (yellow, two end Majorana 
modes, ordering), B (green, one end Majorana mode, ordering), C (red, no end Majorana modes, ordering) and D 
(blue, one end Majorana mode with roles of a and b interchanged compared to phase B, cr^ ordering). The panel at the bottom 
shows the colors of the four phases which would be in increasing shades of darkness in a black and white picture. 

JxJy -Jy = ^^ly-v'^. Using Eq. we get J [cos(20i) + sin(2(/)i) - 1] = M [cos(2(/)2) + sin(202) - 1] For J/M = 1, 

this gives the relationship 

sin(20i +7r/4) = sm{2(f>2 + tt / i) . (34) 

This implies that either (j)2 = ipi+mT or 4)2 + 4>i = n0-|-7r/4, where n is an integer. These describe some of the straight 
line phase boundaries that one sees in the figure for J/M = 1. Some other straight lines in that figure come from the 
conditions that a single root lies at A = ±1 (corresponding to fc = or tt); these give the conditions v + Jx^ IJ. + Jy = 
and V — Jx + ^i — Jy = Q respectively which, for J/M = 1, give sin(0i + 7r/4) ~ ± sin((/)2 + 7r/4). These give some 
additional straight lines corresponding to 4>2 + 4>i — 2n7r ± 7r/2. 
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We have thus provided an exphcit model of long-ranging hopping which has rich spin features and phase diagrams. 
The model supports phases that go beyond those of the Kitaev chain in their multiple Majorana structure. 



V. BROKEN TIME-REVERSAL SYMMETRY: COMPLEX HOPPING AND A 

In this section, we consider a different generalization of our model, one in which TRS is broken. As discussed in 
Sec. II B, this model belongs to a different symmetry class from that of the Kitaev chain, namely class D. We find 
that the fate of the zero energy modes, topology and Majorana structure is unusual in that TRS breaking yields a 
finite regime in the phase diagram that is gapless. 



A. Model and Phases 



Our model is described by the Hamiltonian 

^ = E [ - - ^Vl+i/n + A/n/«+i + A*/I+i/I - Kflfn - 1/2)] , (35) 

n 

where w and A may be complex. We now observe that A can be made real by a phase transformation of /„ and 
namely, changing /„ — >■ /„e'^ and !n^~^^ changes the phase of A by e*^^ without changing the phase of w. 

On the other hand, the phase of w cannot be changed by any transformation without affecting the phase of A. We 
can therefore assume without loss of generality that A is real. Writing w = wqc^'^, where wq is real and positive, we 
obtain 



H = 



J2 [ - Woe^'^flfn+l - Woe-^'f'fl^Jn + Mfnfn+l + fl+jl) - ^i{fVn - 1/2)] . (36) 



By making appropriate phase transformations, we can ensure that (j) lies in the range [0, 7r/2]; we therefore study the 
properties of this model only within this range. If we use the Jordan- Wigner transformation in Eq. p^ . we find that 
Eq. ([5^ takes the following form in terms of a spin-1/2 chain 

H = - [(«;ocos0- A)«+i -f (wocos^-}- A)tT«CT^+i -f- sin0 (ct>J^+i - ct>;^+i) + ^< . (37) 

n 

In momentum space, we find that Eq. (|36p takes the form given in the first equation in Eq. (jl7|) . where 

hk = (2?«osin0sinfc) / — (2i(;o cos^cos/c -f /i) -f (2Asinfc) . (38) 
(Note that we have explicit TRS, h*_f^ — hk, only if (/) = or tt). The energy-momentum dispersion is given by 



Ek = 2wo sin (j) sink ± \J {2wo cos 4> cos k + + 4A2sin^fc. (39) 

The quantum critical points (or lines) are given by the condition Ek = for some value of k lying in the range 
[0, tt], namely, 

(2wo cos (/) cos fc + fj,)'^ + 4A^sin^fc = 4u)o sin^ sin^ fc (40) 

for some k. To understand the complete phase diagram, we look at the quantum critical points for six values of (j) in 
the range [0, 7r/2]. The resultant phase diagrams are shown in Fig. |4l For all values of (except for and tt/2), the 
energy vanishes on two lines and within a two-dimensional region. The two quantum critical lines are given by the 
vertical lines fi/wQ ~ ±2cos(/) (shown in red); on these lines, the energy vanishes at fc = or tt. The two-dimensional 
region (shown in blue) consists of a rectangle whose left and right sides are capped by elliptical regions. The rectangle 
is bounded by the horizontal lines given by A/wo = ±sin0 and the vertical lines fi/wo = ±2cos(/). Finally, the 
elliptical regions on the sides of the rectangle always touch the points {n/wo, A/wq) = (±2, 0). 

So, to describe the evolution from the TRS preserved case {(j) = 0) to the full broken case (0 = 7r/2), at </> = 0, we 
have the phase diagram respected by the Kitaev chain (the first figure in Fig. |4]). Topological phases exist vertically 
below and above the gapless line running from (/^/wg, A/wq) = (—2,0) to (2,0) while the region l/i/wol > 2 are non- 
topological. As seen in Fig. 31 as (j) deviates from zero, the vertical phase boundaries {fi/wo = ±2cos0) approach one 
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another, the gapless hne along A = grows to a two-dimensional gapless regime that consists of a rectangular region 
between the vertical phase boundaries and two elliptical caps on the sides. By invoking adiabaticity, and as shown 
rigorously below, we find that the topological regions persist above and below the gapless regions, and are bound 
between the vertical phase boundaries while the rest of the gapped regime remains non-topological. Finally, upon 
reaching (/> = 7r/2 (the last figure in Fig.|4]), the two vertical lines merge into a single line given by /i = 0, shrinking the 
topological region to non-existence and the two-dimensional region forms a single ellipse given by fj.^ + 4 = Awq . 

The gapless regions, shown in blue in Fig. 21 are quite unusual. For each value of k (except and tt) and (f>, Eq. (j40p 
describes an ellipse in the plane of fi/wo and A/wq- The gapless regions arise when all the ellipses for different values 
of k are combined. This is reminiscent of the Kitaev model of spin-l/2's on a hexagonal lattice; that model also has 
a gapless region in which the energy vanishes at different points in the two-dimensional Brillouin zone [57| . 




<3 




FIG. 4: Phase diagrams of Eq. (|36p as a function of n/wo and A/wo, for (p = 0, tt/IO, it/5, Stt/IO, 2n/5 and tt/2. The 
system is gapless on the vertical red lines and everywhere in the blue shaded regions. The figures show the phases I and II 
which are topological and phase III which is non-topological; the three phases exist in all the figures except for the bottom 
right figure {(p = n/2) where phases I and II do not exist. 



Turning to the fate of zero energy end modes, consider Eq. ([36]) expressed in terms of Majorana fermions: 



H ^ - ^ ^ l^u-ocos^ (a„6„+i + a„6„_i) - A(a„6„+i -a„6„_i) + fi a^bn + u-osinc/) (a„a„+i -I- 6„&„+i) . (41) 



(Based on the time-reversal transformation of a„ and 6„ discussed in Sec. IV B, we again see that that this Hamiltonian 
does not have TRS unless = or tt). The Heisenberg equations of motion following from Eq. (|4T|) give 

wq cos (f) {bn+i + bn-i) - A(6„+i - bn-i) + ^ bn + sin (a„+i - a„_i ) = 0, 
wocos0(a„+i + a„_i) - A (a„_i - a„+i) + ^ a„ + wq sin (/) (6„_i - 6„+i) = 0, (42) 

for the zero energy modes. Eqs. show that the zero energy modes couple a„ and 6„ if sin0 ^ 0, i.e., if TRS 
is broken. In terms of Fig. 01 we find that the the region lying between the two vertical lines but excluding the 
rectangular region forms a topological phase. In this phase, a long chain has one zero energy Majorana mode at each 
end; these modes have real wave functions (involving both a„ and 6„), and they are separated from all the other modes 
by a finite energy gap. We have shown this using numerical calculations but we can also understand it analytically 
as follows. 
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Consider a general quadratic Majorana Hamiltonian of the form 

2N 

H = i CmMmnCn, (43) 

m,n— 1 

where M is a real antisymmetric matrix; hence iM is Hermitian. (For instance, this is the Hamiltonian we get for a 
A^-site system if we define the operators Cji in terms of ct^i and bn a-s C2n—i — s-nd C2n — 

for 7i^l,2,--- ,N). One 

can show that the non-zero eigenvalues of iM come in pairs ±Aj (where Xj > 0), and the corresponding eigenvectors 
are complex conjugates of each other, Xj and cc*; this is because iMxj = Xjxj implies iMx* — —XjX*. The number 
of zero eigenvalues of iM must be even, and one can choose those eigenvectors to be real: if iMxj = 0, we have 
iMx* = 0, and we can then obtain real eigenvectors by taking the linear combinations Xj + x* and i{xj — x*). 

Now, let us consider a Hamiltonian with TRS which has, say, q zero energy Majorana modes at the left end of a 
long chain. Let us assume that the bulk modes are gapped at zero energy, i.e., there are no bulk states with energies 
lying in the range [—Eq,Eo], where Eq is a positive quantity. Next, let us add a small TRS breaking perturbation to 
the Hamiltonian. By adiabaticity, the q zero energy modes at the left end continues to lie within the bulk gap, but 
they need not remain at zero energy. However, we know from the arguments in the previous paragraph, that they 
can only move away from zero energy in pairs. Hence, if q is odd, one Majorana mode must remain at zero energy 
with a real eigenvector. Since we know from Sec. H A that the case with (p = 0, which has TRS, phases I and H 
are topological and have one Majorana mode at each end of a chain, this must continue to remain true if we make 4> 
non-zero, as long as the bulk gap does not close. However, the nature of the end Majorana mode changes from the 
case with TRS (where it involves only a„ or only 6„) to the case with broken TRS (where it involves both a„ and 6„). 



B. General Quadratic Majorana Hamiltonian and Topological Phases 

As seen above, TRS breaking results in couplings of the a fermions amongst themselves, and similarly couplings 
amongst the b fermions. In this section, we will consider the most general Hamiltonian which is quadratic in terms of 
the Majorana fermions a„ and 6„, thus generalizing the arguments presented in Sec. IV A (on long-range hopping) 
to include TRS breaking. This enables us to discuss the Z2-valued invariant which appears in such a system. 

For an infinitely long chain, such a Hamiltonian can be written as 

oo oo oo oo 

= - I ^ ^ Jr ttnbn+r " * X/ X/ l^r anttn+r + Lr bnbn+r], (44) 
r— — oo n— — oo r— 1 n—~oc 

where Jr, Kr, Lr arc all real parameters. In momentum space, this can be written as in the first equation in Eq. ()17p . 
where 

oo 

/ifc = -2 ^ [Jr cos{kr) + JrSm{kr) s"] 

r— — oo 
oo 

2 ^ [{Kr + Lr) sm{kr) I + {Kr - L,.) sm{kr) s^]. (45) 

r=l 

The energy-momentum dispersion is given by 

oo 

Ek = 2 ^ [{Kr + Lr) sm{kr) 



± 2 



Jr cos{kr)\ + i Jrsm{kr)\ + [ ^(A%. - L^) sin(A:r) ] . (46) 



\r— — oo / \r— — oo / \r— 1 



We now observe that Eq. (|45|) has TRS, i.e., /il^. ~ hk for all fc, only if Kr = Lr = for all r. In that case, as 
indicated in Sees. IV A-B, hk defines a vector in the y — z plane as in Eq. (jl9p : since this plane contains the origin, we 
have already seen that we can define a winding number around the origin as in Eq. ([6]), and this is a TI taking values 
in Z (the set of all integers). However, if some of the Kr or Lr are non-zero, then hk does not have TRS and has 
all four components (along J, s^, and s^) in general; one cannot define a winding number around the origin for a 
closed curve (generated by k going from to 2tt) if the curve is in more than two dimensions. On the other hand, at 
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k = and tt, only has a component along given by /i(0) = —2^,, Jr and /i(7r) = — 2^^(— l)''Jr respectively. 
Thus, for the gapped regions in parameter space, we explicitly sec that we can employ the bulk TI defined in Eq. ([7|), 
yielding 



This invariant takes the values —1 and 1 in the topological and non-topological phases respectively; its value cannot 
change unless either h{0) or h{Tr) crosses zero in which case the energy Eq or i?^ must vanish. As an example, 
for the simple system with complex hopping discussed in the previous subsection, Sec. IV Al we see that Vbuik = 
sgn(/i^ — Awq cos^ 4>) is equal to —1 in the regions \fi\ < 2wo cos (f>. Once we consider the fact that in this regime, the 
system is gapped only for |A| > 2wosin0, the condition Vbuik = ^1 captures the correct topological regions shown in 



To explicitly analyze the nature of the end Majorana modes, we have seen that in a system with TRS (with 
Kr = Lr = 0) there can be any number of such modes at each end. For instance. Fig. [1] (e) shows a system in which 
only J2 is non-zero and there arc two zero energy modes at each end of an open chain shown by 03, 04, 65 and fog. 
However, if we now break TRS slightly by turning on small values of Ki and Li in the Hamiltonian, i.e., introduce 
terms like —iKia^a^ — iLib^be, the energies of these modes shift away from zero to ±ifi and ±Li, destroying any 
zero energy Majorana modes. One can see an odd-even effect manifest in that if there are an odd number of zero 
energy Majorana modes at one end of the chain in the presence of TRS, breaking the symmetry would in general 
couple these modes. But given that there are an odd number of states and that the Hamiltonian only permits states 
to come in ±E pairs, unlike in the even case, one E = Majorana end mode has to survive. 

To conclude, a system with TRS can have an arbitrary number of zero energy Majorana modes at each end of a 
long chain. If TRS is then broken, pairs of these end modes move away from zero energy. In general, therefore, we 
are left with either no Majorana modes or only one Majorana mode depending on whether the parent system with 
TRS has an even or an odd number of Majorana modes. 



Having considered the effects of long-range hopping and TRS breaking on the topological features of the Kitaev 
chain, we now treat the effect of subjecting the chain to a spatially varying potential landscape. For the discretized 
lattice, the potential landscape can be described in terms of a site-dependent chemical potential. The last term 
containing the global chemical potential ^ in Eq. ([1]) can be replaced by one having a local chemical potential 
The topology of such systems has been studied before [H, [H, M, H H M, IHMI, Hi 

The phase diagram in Fig. [5] (a) shows that in the the spatially homogeneous case, topological phases exist in the 
Kitaev chain even with arbitrarily small superconductivity present (A ^ 0) and persist until a chemical potential fi 
that is of the order of the bandwidth w is applied [l|. Here we shall see that this is a singular limit and that for 
spatially varying potential these features no longer persist [36l. [stI. [50| . In fact, in general wc will see that for any finite 
potential, a minimum amount of superconductivity is required to reach the topological phase, though there remain 
interesting counterexamples to this general rule (see Fig. [TU)) [s^ . 

In what follows, we analyze three different kinds of potential landscapes, expanding and building on our recent 
short work, Ref. (30|. As the simplest example and a direct application of the transfer matrix techniques presented 
in Sec. II B, we obtain the phase diagram for periodic potentials having short periods commensurate with the lattice 
(see also [U, IH, (H). We then develop our transfer matrix further to enable an extensive study of more complex 
potentials. In particular, wc provide a mapping between the normal state properties of the wire and its topology in 
the superconducting state. Our result is a more general form of the mappings presented in Refs. [27l. [28|. However, 
this transformation becomes singular at A = 1, and thus cannot by itself be used to obtain the full phase diagram 
(see Sec. VI. B). The line A = 1 may be treated separately (this regime corresponds to the random field Ising model). 
We have also developed a mapping between the topology between the gap function A and 1/A. Taken together, 
these procedures allows us to obtain the topological phase diagram and the decay length of the end Majorana wave 
functions in terms of the normal state localization properties of the system. This mapping enables us to leverage the 
vast literature on the localization properties of ID systems to the question of the topology of these systems in the 
presence of superconductivity. 

Given the close connection between normal state properties and topology, we find that a great deal of insight can 
be gained by briefly reviewing the band structure of simple periodic systems, which is done in Sec. VI A. In Sec. 
VI B, the mapping between normal state properties and topology will be developed. Finally, Sec. VI C presents an 
application of these methods to ultra-long periodic and quasiperiodic potentials. 




(47) 



Fig.m 



VI. SPATIALLY VARYING POTENTIALS 
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FIG. 5: Topological phase diagram showing the topological (T) and non-topological [N) phases as a function of the potential 
strength W and the superconducting gap A for (a) a uniform potential /i„ = W, and (b) a periodic potential having the pattern 
iW,W,-W,-W). 



A. Periodic Potentials and band structure 



Given that we will forge a connection between the topological and normal state prop erties of inhomogeneous systems, 
it is helpful to review some basic properties of the normal state band structure [30l . [sij . The problem of an electron 
hopping in a periodic potential is described by Eq. ([1]) where we set A = and take the chemical potential /i„ to 
depend on n in a periodic way. For a state with energy we get the discrete form of the Schrodinger equation 

- w(/„+i + fn-l) - fJ.nfn = Efn, (48) 

This can be written in the transfer matrix form 




1 




(49) 



The energy spectrum consists of those values of E which admit plane wave states, /„ ^ giifen^ p^j, system with 
period q, the appropriate object to consider is 




(50) 



Since ^ is 2 x 2, its eigenvalues satisfy 



- (TiA) A + dctA = 0. (51) 
Let A± denote the two eigenvalues of A. Since det ^ = 1, we have A+A_ = 1. Further, A+ -|- A_ must be real. Hence 

'e±*'' for |Tr A\ < 2, 
A± = { e±'' for Tr A > 2, 

e±'' for Tr A < - 2, 

where r is real. The object Tr ^ is a g*'* order polynomial in E, and the spectrum of the system corresponds to the 
those values of E for which |Tr^| < 2. A useful quantity which quantifies the localization properties of the system 
is the Lyapunov exponent j{W) = \imj^^^ jj In | max(A-|-)|. Fig. 5 presents a detailed example of the normal state 
properties of a periodic system at zero energy as well as anticipating the connection between the normal state and 
the topological properties developed in See. IV B. 



It is instructive to consider the case of simple periodic patterns as they offer a great deal of insight into the properties 
of the phase diagrams in more complicated situations. These features anticipate and provide qualitative understanding 
of our central results in subsequent sections. For instance, we will consider potentials of the form = Wti, where 
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TABLE II: Criteria for topological phases for a selection of periodic potentials (we have set w — 1). 



period 


pattern 


topological for 


1 


...,w,w,w,... 


\W\ < 2 


2 


...,w,-w,... 


A > \W\/2 


4 


...,w,w,w, -w, . . . 


> - 1 


4 


...,w,w, -w, -w, . . . 


A > WV4 



is some inhomogeneous potential of fixed strength and zero mean. A (very) naive caricature of such a system would be 
a period-2 potential of the form /i„ = —W, W, ... or a period-4 potential fin = —W, —W, W, W, .... The topology 
of these simple periodic patterns may be obtained by taking the transfer matrix A = Jln^i •^n introduced in Sec. II 
B, where q is the period of the potential. Given A, TI i> may be straightforwardly found by applying the methods of 
Sec. Ill B~2] and Eq. (|12p. in particular. Such simple periodic patterns naturally arise in a spin ladder model in which 
a vortex degree of freedom controls the sign of the chemical potential in the corresponding fermion system (sij . a 
one-dimensional version of Kitaev's celebrated honeycomb model [STf . Several simple examples arc given in Fig. [5] 
and Tabic ini)- 

The presence of a phase offset in commensurate periodic potentials can significantly affect the topological phase 
diagram. As an illustrative example, let us consider the period-4 pattern given by 

fin = V2ficos(^^{2n + l) + . (52) 

Since shifting cj) ^ (j) + tt/2 is equivalent to shifting n n + 1, it is sufficient to study this problem in the range 
< < tt/2. Upon finding the eigenvalues of the product of four transfer matrices A4A3A2A1 of the form given in 
Eq. ((54)) below, we find that the phase boundary between the topological and non-topological phases is given by 

4V2|wA| = /iVl + cos(4(/)). (53) 

The phase diagram is shown in Fig. |6l The fact that the boundary between the topological and non-topological phases 
diverges in as — > 7r/4 can be understood as follows. For = 0, we know that the pattern of the chemical potential 

takes the form P = + H ... and thus has the phase boundary given by /i = 2-\/|?iiA| (see the last line of Table 

I). For (f> = 7r/4, the pattern assumes the form Q = +0 — + .... In the limit — > 00 (in which the hopping plays 
no role), the eigenstates of the system are localized to each site. Given that ^ = at every other site, the system 
must contain bound Majorana modes and therefore must be topological. The general shape of the phase boundary 
then follows from the fact that it connects the points (2-\/|?iiA|, 0) and (oo,7r/4) in the {fJ,,(f>) plane. In general we 
conjecture that for any periodic pattern, fin = fi co&{'KQn -\- (p) where Q is a, rational number, if the chemical potential 
vanishes at certain sites for some value of 0, a 'finger' of the topological phase will extend up to = 00 for that value 
of (t>. As will be seen in the next section, the sensitivity of the phase diagram to potentials vanishing at certain sites 
is quite ubiquitous and results in a non-analytic behavior of the phase boundary for disordered systems. 



B. Mapping to Normal Systems and Duality 



As observed in Ref. [5Jl, the product of such transfer matrices, A = strongly reminiscent of that used 



to determine localization properties of the normal state Anderson disorder problem, i.e. Eq. (|49|) . Wc build on this 
observation by explicitly mapping the topological phase diagram to the normal state properties of the system (similar 
to a procedure developed in |28|). This mapping is performed in two steps. The first allows us to determine topological 
regions in the superconducting system in the range < A < 1 based on a knowledge of its normal counterpart (same 
spatially varying potential but no superconductivity). The second is a duality which relates the topological phase 
diagram for a given superconducting gap A to that corresponding to strength 1/A. Taken together, these provide a 
complete mapping between the normal state and the topology of the superconducting system. 

Mapping to Normal Systems - In Sec. II B, wc introduced the transfer matrix relevant to a superconducting 
wire having spatially varying potentials: 

/ Un A — lU \ 



An = ( ^+^" ^J'" 1 . (54) 
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FIG. 6: Topological phase diagram for a period-4 potential, showing the topological (T) and non-topological (A'') phases as a 
function of fi and (j). 



For < A < 1, we perform a similarity transformation An ~ VSSAnS ^ with S = diag((5^/'*, and S = 

The matrices A„ are of the form shown in Eq. ([TU]) with A — >■ and /i /i„/-\/l — A^. This immediately gives 




A{W, A) = 



Taking the logarithm of the eigenvalues of Eq. ((55)) . the condition that IA2I = 1 is given by 



W 



r,0 



1- A 



(55) 



(56) 



reminiscent of a result in [2^, where we have defined the Lyapunov exponent 7(W,A) = \imj^^^ In |A2(W^, A)|; 
the Lyapunov exponent is the inverse of the localization length, 7(W^, A) = l/i(W,A). In the limit 7(IV, A) 0, 
the system is gapless and Eq. ([55)1 describes the phase boundary separatin g th e topological and non-topological 
regions of the phase diagram. This relation quantifies the observation in [36l. Ts^] that in general a critical amount 
of superconductivity must be applied before the system is driven into a topological phase. For the case in which the 
system is metallic (i.e., 7(VF, 0) = 0), any non-zero A will give rise to a topological phase. 

If /i„ = for all n, we see from Eq. ([M]) that all the An are identical and have both eigenvalues less than 1 in 
magnitude if A > (we are assuming that w > 0). This implies that ^{W = 0, A) < 0. Thus the line W = and 
A > will always lie in the topological phase. 

Duality - The similarity transformation (Eq. (|55p ) is only valid in the range < A < 1. The form of the phase 
diagram for A > 1 may be obtained by noting that the transformation 



where P — ; 



A 1/A, and P 
P: ^ {(-1)>4, 



(57) 



leaves the eigenvalues of A unchanged for Af even. Thus, if a point (Wo, Aq < 1) lies on the phase boundary of P, then 
(Wo/Aq, I/Aq) lies on the phase boundary of P. This duality strongly constrains the form of the phase boundary in 
the cases in which the distribution is invariant under the transformation in Eq. (j57p . 
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FIG. 7: (a) Plot of the Lyapunov exponent for a period 10 potential with Hn — /j, for n = mod 10 and — /i otherwise. 
Inset shows Tr A. The allowed plane wave states (corresponding to 7?(/i) = 0) have |Tr A\ < 2. (b) The ^-A topological phase 
diagram (topological region indicated in white, non-topological region in gray). Note that the values of ^ for which 7(/i) = in 
(a) are precisely those regions which become topological for arbitrarily small A. This relationship is elucidated in Sec. VI B. 

An interesting illustration of this duality is provided by the periodic patterns in Table [III For example, note that 
the uniform case (period 1) and the period 2 case are dual to each other. Indeed, taking A — >■ 1/A and W — > W/A 
for the equation W = 2 (phase boundary for period 1 case), we obtain the phase boundary A = VF/2, appropriate for 
the period 2 case. Similarly, the two period 4 patterns are self-dual. For instance, taking A^ = — 1, and applying 
A ^ 1/A and W -> W/A gives l/A^ = W^/A'^ - 1 which is equivalent. The same holds true for A = W^/A. 

Random- field quantum Ising Chain - Finally, at the point A = 1, the system maps to the quantum Ising 
chain subject to a spatially varying transverse field. As can be seen in Sec. Ill C[ the Jordan- Winger transformation 
applies even when the local chemical potential fin has a spatial dependence and appears as a spatially varying field 
along the ^-direction. Along the line A/w = 1, we have Jy = 0, and thus for a random distribution of fin, the system 
corresponds to the ID random transverse field Ising model, which has been studied in great depth [2y, [H^. The 
matrix A{W, 1) has the eigenvalues ^ O^i ^'^'^ 0. Eq. ([T^ reveals that the phase boundary passes through the 
point for which 

(ln|/i„|) =ln2, (58) 

where (In |yu„|) = limA^^oo jf J2n=i In |Mn|• 



C. Ultra-long Period and Quasiperiodic Potentials 



Building on our work in Sec. VI A, we consider periodic potentials of the form /.j„ = B cos (27rcjn — 9) but now 
allow the period of the potential (g, for uj = 2Tip/q, for p,q £ N)) to be either very long or infinite p^ . [35| . The 
Hamiltonian formed from this potential is known as the almost Mathieu operator. For B = 2, the resultant equations 
of motion are known as Harper's equations and arise in the context of a charged particle moving on a 2D square 
lattice in the presence of an external magnetic field [soj . 

For UJ irrational, the associated normal state energy spectrum is a fractal known as Hofstadter's butterfiy [59j . 
This fractal structure arises from a property of the quantity Tr A (with A = 0, see Sec. VI A) for the operator 
fin = 2cos(27rwn — 9) (see Eq. ([TU])). This operator has the special property that for uj = p/q,p,q E N, there are 
q distinct bands [5^. This is a special property of Harper's equation. For instance, note that in Fig. [3 although 
q = 10, only 5 distinct bands exist. In contrast, Fig. |S] shows that all g = 10 bands are present for Harper's equation. 
For example, if we imagine letting uj — >■ \/2 by taking a sequence of better and better rational approximations, the 
number of bands will grow without bound. This is the origin of the unusual point-set topology of the spectrum for 
irrational u. A plot of the allowed bands as a function of the energy and uj leads to a fractal known as Hofstadter's 
butterfly [s^. The general features of Hofstadter's butterfly maybe seen in Fig. EJa). Indeed, for uj irrational, the 
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FIG. 8: (a) The Lyapunov exponent 7 of the normal state (A — 0) for fin = V + 2 cos {2TTLon) with to = 1/10. (b) Topological 
phase boundary showing the merging of the topological regions as described in the text. For A ^ 1, there are are 10 distinct 
regions which are topological. At A = 0.2, the four central regions have merged to form a single region. 

total 'length' associated with the spectrum is zero (by length, we mean J f{V)dV, where / is 1 if F corresponds to 
a point in the spectrum, if not). Formally, the Lcbesgue measure of the almost Mathieu operator is |4 — 2|i3|| [60| . 
Very roughly, as g -> od for B ~ 2, the length of the spectrum (as a function of energy) is zero, although there are still 
an infinite number of distinct bands! These seemingly paradoxical properties are the features of a so-called Cantor 
set. 

These remarkable properties lead us to consider the topology of a ID system subject to a potential of the form 

fin + 2cos{2npn/q) . (59) 

The allowed zero energy plane wave states in the uj V plane correspond to Hofstadtcr's butterfly. Of interest are 
the properties of the topological phase diagram as A is 'turned on'. This problem is easily solved given the Lyapunov 
exponent 7 of Harper's equation which we have plotted in Fig. IHJa). This plot is different from the usual plot of 
Hofstadtcr's butterfly in that it shows the value of the Lyapunov exponent 7 rather than just the spectrum (the 
spectrum corresponds to points for which 7 = 0). Indeed, a feature of this plot are characteristic horizontal striations, 
which show the sensitivity of the localization length to the period of the potential (g) . 

We begin by considering the topological properties of a potential with a given u. Fig.|8]plots the Lyapunov exponent 
of a particular potential {lu = 1/10) as a function of V. As shown in Fig.[51[b), for the same potential and A ^ 1, there 
are q topological regions inherited from the normal state. These distinct regions fuse as A is increased. The precise 
value of A for which two topological regions merge is determined by the strength of the Lyapunov exponent between 
the gaps. Since this quantity tends to be larger for larger \V\, the gaps closer to ^ = tend to merge before those 
with larger \V\. We now generalize to the full ut — V space. As expected from our general analysis, the normal state 
properties (Fig. [9lja) directly inform the topological phase diagram (Fig. [9jb). We see that as superconductivity is 
increased, Hofstadtcr's butterfly is 'filled in' by non-topological regions. The value of A at which the system becomes 
topological is extremely sensitive to the period q of the potential. 



The case of Harper's equation allows us to study the crossover exhibited from periodic to quasiperiodic potentials. 
We now consider the case Hn = W cos (27ra;n) with oj being irrational. The normal state features of this potential 
have been well-studied; the system is metallic (i.e., there are plane wave states at = 0) for W < 2 and exhibits 
a metal-insulator transition at the critical value W 



-yiW, 0) = In {W/2) for W 
phase for 



> 2 and for < < 2 



2 

for w 



61| . The normal state Lyapunov exponent takes the form 
irrational [g^, [11] . Eq. ([55]) then predicts a topological 



A> -W 
2 



1. 



(60) 



This result holds for all values of A > given that the transformation uj — )■ 
duality transformation, A — > 1/A and W — >■ W/A, leaves Eq. ([50)) invariant, 
point (W, A) = (4, 1) lies on the phase boundary. 



a; -I- 1/2 yields Eq. and that the 
Finally, Eq. ([55| also shows that the 
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FIG. 9: (a) A colorscale plot of 7(V^, 0) for u — n/200, < n < 200. Darker regions correspond to smaller values of 7. The 
characteristic striations show the spectrum's sensitivity to the period of the potential; i.e., small changes in uj can lead to large 
changes in the period (which is given by q, where lj = p/q, with p and q relatively prime). The resulting figure is reminiscent 
of the fractal known as Hofstadter's butterfly, (b) The topological phase diagram for A = 1/5 mimics the low- lying values of 
7 in (a). 




FIG. 10: Topological phase diagram for a potential Hn 
non-topological regions respectively. 



W cos {2TTLun) . The symbols T and A'' refer to topological and 



VII. DISORDERED POTENTIALS 



The topic of disordered superconducting wires that break the SU(2) symmetry associated with spin, namely those of 
the symmetry classes D and BDI discussed in previous sections, has been actively researched for over a decade (47ll49l- 
[5]| . One of the highlighting features of these systems is that their symmetry properties greatly alter localization 
physics. In particular, while Anderson localization dictates that states are always localized in one-dimension in the 
presence of even the weakest disorder, these systems allow for the presence of a critical disorder point in parameter 
space that permits a delocalized state at zero energy [U ligj - fsij . The point acts as a mobility edge in that it separates 
two localized regimes. In these systems, a variety of approaches have probed the manner in which the density of 
states diverges at zero energy, odd-even effects for coupled chains, and conduction properties (since charge is not a 
conserved quantity for excitations about the superconducting ground state, one studies thermal insulating/conducting 
properties). In pioneering work by Motrunich et al. [s^l, it was shown that the localized phases separated by the 
delocalized point are fundamentally different from one another. One phase supports end Majorana modes while the 
other does not, consistent with the classification of the phases as topological and non-topological, respectively, as in 
the disorder-free case. 

Recent studies have focused on this feature of the presence of Majorana edge modes in the context of disordered 
systems 0, H^-IH, US HE 113, Along these lines, here we exploit the transfer matrix techniques that we have 
developed in previous subsections to pinpoint the conditions for the existence of these end modes in a full range 
of models of disorder and disorder strength. The map made to normal systems in Sec. VI B allows us to borrow 
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FIG. 11: Topological phase diagrams for (a) box and (b) double box potentials. The red dotted lines indicate that at A = 1, 
the slope of the phase boundary, dH^/dW , has a discontinuity in (a) but is continuous in (b). 

extensively from literature on Anderson localization and perform a comprehensive study. We derive general features 
of the phase boundary for disordered superconducting systems, including a characteristic discontinuity of the phase 
boundary at the point associated with the random field Ising model (A — 1). We also obtain topological phase 
diagrams for a variety of disorder distributions based on our mapping to an extensive range of known results from 
Anderson localization studies of normal systems. 




A. Setup and General Results 

We consider Eq. ([T]) with replaced by a spatially dependent The /Xn is drawn from some distribution and is 
uncorrelated, i.e. {unfJ-n') = USn^n' and (/i„) = 0. The quantity U is the standard deviation of the disorder and thus 
characterizes its strength. 

For weak disorder, the Lyapunov exponent may be obtained from perturbation theory for the normal state sys- 
tem [sl, [fil] and is given by 



2 



Applying Eq. ([5^ . we obtain the condition for the topological phase 



A > 



This result may be compared to that of a continuum model based on the Dirac equation, which gives a non-topological 
phase for A > = 0.125J7 (see Ref. ^^). For disorder distributions that are symmetric around 0, the self-duality 

condition P = P in Eq. (|57p is satisfied. In this case, we can employ this duality transformation, A — >■ 1/A and 
U U/A'^, to show that Eq. (p^ also describes the phase boundary in the limit of strong disorder. Finally, we 
mention that near the Ising point (A = 1), the disorder may give rise to a discontinuity in the phase diagram. We 
describe the physics of this phenomenon in the examples given now. 



B. 'Box' Disorder 



As a generic representative for disorder, we now turn to the case of 'box' disorder for which the probability of 
fin falling at any point in the range —W/2 < < W/2 is equally likely. The low-energy behavior as shown in the 
numerical simulation in Fig. fTlT a) is in good agreement with Eq. (for box disorder, U = W'^/12). Eq. ([55)1 reveals 
that the phase boundary passes through the point {W, A) = {Wc, 1), where Wc = 4e « 10.873 (box disorder) with e 
being the base of the natural logarithm. 

Discontinuity of phase boundary- A noteworthy find is the observed discontinuity suffered by the phase 
boundary as it passes through the random transverse Ising field point A/w = 1 discussed in the previous section 
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(Figs. [TTJa) and fT2)) . We can calculate the phase boundary near A/w w 1 by noting that according to Sec. IVIBl 
^ ^ ^^^2 , 0^ , the effective Lyapunov exponent that we seek in Eq. (j56p . corresponds to that of very strong disorder 
for A — > 1. In this limit, we can use the known form of the normal state Lyapunov exponent for W ^ 1 

-/{W) - In (w/2 + ^1^2/4 -4) - ^M^2/4-4 - In 2, 



In {W/2) - 1 + A/W^ + O (l/W^) 



(63) 



(there is a typo in the expression given in (65|). Substituting this expression into Eq. ()56p and invoking self-duality, 
we obtain the phase boundary (to linear order around {W, A) — (4e, 1))) 



A 



2e2+2 



2e2 



w- 



4 for A < 1, 



4^ for A > 1. 



(64) 



As seen in Fig. Illf a). this result is in reasonable agreement with numerical simulations. 

In order to understand the origin of this discontinuity, we investigate the phase boundary by performing perturbation 
theory in the quantity (1 — A^). A straightforward calculation shows that the eigenvalues of the matrix A are and 



1 



N 



(1 + A) 



N 



n ^« 



1 



1 



(65) 



From this expression, it is straightforward to obtain the Lyapunov exponent in terms of the disorder distribution. 
Letting A = 1 + J, we get 



7 =(lnAi)-ln2- 1(1 + 4(1/^)2)5 

up to order . Let us define 8 = we then get that the phase boundary for 8 

(In |/i„|) =ln2- (1 + 4(1^)2)^ + 0(52). 



obeys 



(66) 



(67) 



This shows that the phase boundary is fragile towards singularities when /x„ is allowed to come arbitrarily close to 
zero. Indeed, our simulations have shown that the discontinuity is absent for disorder distributions that avoid zero 
energy. A large class of disorder distributions that cover zero energy give rise to a discontinuity in the slope of the 
topological phase boundary at A = 1. This is another manifestation of the sensitivity that the phase diagram shows 
to /i„'s which are equal to zero (for example, see Eq. (|58p or Fig. 7). 



C. Other forms of disorder 



'Double Box' Disorder- A useful test of our hypothesis on the fragility of the phase boundary for arbitrarily 
small values of /i„ is to examine a system for which the disorder distribution excludes \in = and thus we expect that 
the phase boundary be continuous near A = 1. Consider the following distribution for the local chemical potential, 

(68) 

I otherwise. 

Here, Wc = 2e « 5.437. Note that all the {1/ ji^) exist. We can find the behavior of the phase diagram near A = 1 
using Eq. (|67p . More directly however, we note that to linear order near A = 1, the only expression which is invariant 
under Eq. dST]) is 

A = aVK-l. (69) 

This expression well describes the phase boundary for a = 1/e. This phase diagram is shown in Fig. Illf b) and is in 
good agreement with this prediction, in particular, being devoid of the singularity at A = 1. 
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FIG. 12: Topological phase diagram for Lorentzian disorder. 



Lorentzian Disorder- Finally, we turn to the specific case of disorder which is unbounded and has a diverging 
standard deviation, the Lorentzian case. The distribution of local chemical potential is drawn from a distribution of 
the form 



1 W 
f^{x;W) = - 



The phase diagram is exactly soluble in this case since the normal state density of states is known exactly [& 
zero-energy Lyapunov exponent, first obtained by Thouless [63, takes the form 7l(W^, 0) = In (w/2 + ^/l 
Once again invoking Eq. (|56l) and self-duality of the phase diagram yields a phase boundary 



(70) 



The 



W = 



2A for A < 1, 
2 for A > 1. 



(71) 



This result, as shown in Fig.[T21 is in excellent agreement with numerical simulations. It should be pointed out that the 
features of this phase diagram arc extremely unusual. For instance, Eq. (j62p fails to hold because the second moment 
(/i^) is ill-defined for //,. This example is noteworthy since, for W > 2 the system is always in a non-topological phase; 
no amount of A can drive the system topological. Studying these examples has shown us, among other features, that 
typically the larger the disorder, the more superconductivity is required for Majorana end modes to exist, and that 
the topological phase diagram is highly sensitive to the nature of the disorder distribution. 



We also mention that the method presented here may be applied to a Guassian distribution of disorder where a 
closed form of the Lyapunov exponent is known (see (27j for an analysis of this example). However, the result involves 
Airy functions and is quite complex and we omit this example given that it has the same qualitative behavior as the 
phase boundary for 'box' disorder. 

Future Directions - A very natural next step would be to consider spatially inhomogeneous potentials for a 
system in class BDL The interplay between such potentials and phases exhibiting multiple Majorana modes should 
be quite rich. We offer a hint of this in Fig. [131 The figure indicates the intricate behavior which may occur in such 
cases, including re-entrant phases. Given this exciting behavior, we note that applied potentials offer a very promising 
route for the engineering of topological phase diagrams. We note however that the simple numerical methods which 
were applied here for a disordered system {L ^ 10'*) with nearest neighbor hopping will not work for systems with 
longer range hopping. This is due to the difficulty in extracting all the eigenvalues of a large product of matrices (and 
not just the largest one). 



VIII. SUMMARY AND DISCUSSION 



In this work we have presented a comprehensive study towards enhancing our understanding of topological phases 
of one-dimensional superconducting systems, TIs and Majorana end modes. We began with reviewing a prototypical 
model, the Kitaev chain, which describes spinless electrons hopping with an amplitude w between the sites of a 
lattice, p-wave superconducting order (denoted by a parameter A), and an on-site chemical potential fi. This model 
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FIG. 13: (a) The topological phase diagram of a uniform class BDl system (for which A is constrained to be real) described 
by Eq. (|16|) with Jo — ^/2, Ji = J2 = {w — A)/2, and J_i — J-2 = {w + A)/2, with w = 1. That is, the system is similar 
to Eq. Ul but with next-nearest neighbor hopping and superconductivity equal to w and A. The different regions are labeled 
according to the topologically protected quantity n/, which reflects the Majorana end mode structure as defined in Eq. (|lip . 
The phase diagram was obtained using the method presented in Sec. Ill B 21 (b) The same phase diagram for a 15-site system 
with a yUn = nVn, where Vn is a. random real variable between —1 and 1. 



has attracted a great deal of attention as the simplest model which exhibits topological phases and non-topological 
phases depending on the parameters appearing in the Hamiltonian. In the topological phases, a long chain with open 
boundary conditions has one zero energy mode localized at each end. The non-topological phase has no end modes. 
The model can be mapped, using the Jordan- Wigner transformation, to a spin-1/2 chain with XY couplings and a 
magnetic field applied along the z direction. There are two topological phases which have long-range order in the x 
or y components of the spin, while the non-topological phase does not have any long-range order. 

Next, we showed that there are a number of TIs which can be used to characterize the different phases for uniform 
systems. There are different invariants depending on whether the system is time-reversal symmetric or not, commonly 
referred to as belonging to symmetry classes BDI and D, respectively. The TIs can be defined by invoking the bulk 
or the boundary properties for a large system with periodic and open boundary conditions respectively. For time- 
reversal symmetric systems, the bulk invariant is a Z-valued winding number obtained from a closed curve lying in a 
two-dimensional plane, formed by the set of Hamiltonians spanning momentum space. If time-reversal symmetry is 
broken, the closed curve no longer lies in a plane and a winding number cannot be defined. However, the Hamiltonians 
at the points k = and tt can be used to define a Z2- valued invariant which can take values ±1. As for boundary 
invariants, based on the Heisenberg equations of motion for open chains, we obtained TIs that count the number of 
Majorana modes (whose components we explicitly argued to be real) at the ends of the system as well as their parity. 
We showed that in the presence of time-reversal symmetry breaking terms, multiple Majorana modes, while perhaps 
confined to the ends, disappear by obtaining complex components and moving away from zero energy, thus allowing 
at most a single Majorana mode at each end. We illustrated the bulk-boundary correspondence which relates the 
values of the bulk and boundary invariants to each other. 

While the Hamiltonian for a Kitacv chain involves only nearest-neighbor hoppings and an on-site potential term, 
both only have interconnections between two types of Majorana fermions, a and 6, we generalized this model in several 
ways. First, we introduced long-range hoppings of the Majorana fermions to construct time-reversal symmetric models 
which can have any value of the Z-valued TI in the bulk, and therefore any number of Majorana end modes. For a 
time-reversal symmetric system, we found that each end mode is purely of type a or type b. As an illustration of the 
kind of rich phase diagram that one may have with long-range hoppings, we considered a particular model in detail. 
This model has four terms in the Hamiltonian one of which describes a next-nearest neighbor hopping. The model 
has four phases with the number of end modes of types a and 6 varying from zero to 2. The model can be mapped to a 
spin-1/2 chain; the corresponding model enjoys a complete duality. All the four phases have long-range order in terms 
of X and y components of either the original spins or the dual spins; the phases are related pair-wise to each other 
through duality. Second, we introduced time-reversal breaking hoppings of Majorana fermions on nearest-neighbor 
sites. We again obtained topological and non-topological phases depending on the parameters. The topological phases 
have zero energy Majorana end modes but now each end mode involves both a and b operators. Interestingly, we 



28 



found extended non-topological regions in parameter space having a vanishing bulk gap. It is unusual to have such 
gapless phases in parameter space; the other example of this that we are aware of is the Kitaev model on the hexagonal 
lattice in two dimensions (57| . 

Next, we explored a number of models in which the chemical potential /i„ varies as a function of the site label n. We 
considered three different situations in increasing order of complexity resulting in a rich variety of phase diagrams. In 
order to find these phase diagrams, we developed an important connection between the transfer matrix of models with 
superconductivity and the transfer matrix of non-superconducting models. Namely, we showed that the localization 
length of the end modes of a model with A 7^ can be related to the Lyapunov exponent of a model with A = 0. This 
connection enabled us to make full use of the vast amount of literature available for non-superconducting systems 
with complicated patterns of chemical potentials to identify topological regimes when superconductivity is present. 

The three kinds of patterns of that we analyzed are as follows. First, we studied the case where is a periodic 
function of n commensurate with the lattice. We found that the boundary between the topological and non-topological 
phases depends on the period, amplitude and phase of /j,„. We noted, interestingly, that the boundary displays a 
cusp-like structure if jin happens to vanish at certain sites in a periodic way. Second, we looked at the case in which 
/in is quasiperiodic. Here we discovered a phase diagram that has a fractal structure reminiscent of Hoftsadter's 
butterfly (which appears in the energy spectrum of electrons moving on a square lattice in the presence of a magnetic 
flux through each square which is an irrational multiple of the flux quantum). Third, we explored the case that 
is disordered and is drawn from some probability distributions. Depending on the distribution and the value of the 
superconducting order parameter A, we find a wide variety of phase diagrams. If the disorder distribution includes a 
range of values around /!„ = 0, we found that the phase boundary has a discontinuity in the slope at the point where 
A is equal to the hopping amplitude w. We showed this analytically using perturbation theory in the parameter 
A/w — 1 and we confirmed this numerically for the box and Lorentzian distributions. In contrast, there is no such 
discontinuity at A = if the distribution does not allow the possibility of /i„ = 0, as we showed for the double 
box case. (In order to facilitate these studies, we developed a duality relation between regions with A/w > 1 and 
< 1). The relation between the superconducting system and the problem of Anderson localization in normal electronic 
systems enabled us to analytically find the topological phase diagram in some cases (the quasiperiodic pattern and 
Lorentzian disorder), and to numerically find the phase diagram in other cases in an efficient manner. 

In the future we can consider many other generalizations of our work. As mentioned at the end of Sec. VI, a 
very natural next step is to consider multiple Majorana modes in disordered systems obeying time-reversal symmetry. 
Additionally, we can study systems with spinful electrons [47, 49 51, j)^] and multi-channel wires [6 IjJ. Interactions 
between electrons are known to play a key role in the behavior of one-dimensional systems, and it would be useful 
to study the effect of interactions on the Majorana end modes [1, [l^ - [l7| in the generalized systems presented here. 
Finally, the effects of finite temperature (through electron-phonon interactions) and finite system length (which will 
lead to mixing between the end modes) also need to be studied in detail. 

In part, the interest in realizing Majorana fermions in a condensed matter set ting is motivated by the prospect 
of performing topological quantum computation by manipulating these particles (68[. Recent theoretical work has 
shown (4^ that the manipulation of Majorana modes in a T-junction geometry enacts unitary transformations on the 
(nearly) degenerate subspacc formed from the Majoranas and the Dirac fermions they share. A crucial clement of this 
proposal is the need to dynamically control the topology of segments of the wire. In the case of a clean uniform wire 
(as was seen in Sec. 11), the magnitude of the chemical potential alone controls the topology. Of particular relevance 
to such proposals, we have shown that applied potentials allow the superconducting gap to serve as another 'knob' 
for the topology, and thus an engineered potential could be used to exert more precise control over the topology of 
segments of the system. Many of the theoretical tools introduced here may be applied to questions of the physics deep 
within the topological phase, not just at the phase boundary. Topological quantum computation requires a careful 
consideration of the coupling between different Majorana bound states, crucial for the initialization and read-out 
steps in the protocol as well as understanding the effects of unwanted evolution between states in the ground state 
subspace [6l|- The analysis carried out in Sec. VI to determine the topological phase boundary may be easily modified 
to find the localization length of the Majorana mode anywhere in the topological phase. This provides a means of 
estimating the coupling between Majoranas. 

Finally, we turn to the possibility of testing our results in experimental systems. The effectively spinless models 
that we have studied in this paper are generalizations of the experimental proposals of Refs. in the limit in which 
the Zeeman energy greatly exceeds the superconducting gap and the spin-orbit energy scale, and are applicable to 
the setting of the pioneering experiment of Ref. [11] . This experimental setting thus provides a natural playground to 
explore our generalizations, particularly in the the context of spatially varying potentials. In principle, the predicted 
rich slew of topological phase diagrams can be explored by a controlled application of periodic and disordered poten- 
tials. The effects of time-reversal symmetry breaking can also be studied in a controlled fashion by the application of 
magnetic fields A highly exciting prospect would be to discover phases with multiple independent Majorana modes 
at the ends of a wire. While previous studies showed the existence of such phases in multi-channel wires, perhaps 
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potentially realizable by fabricating a system involving several coupled wires, we explicitly show here that they could 
exist within a single wire in the presence of long-range hopping. 

The detection and study of exotic particles was once thought to be the exclusive purview of high energy physics. 
If the early experimental results are borne out, Majorana fermions give us yet another example of the synergy 
between quantum mechanics and many-body effects, rescrambling electrons whose identity is so commonly thought 
inviolable, into new and exotic objects. Our study of the numerous ways in which Majorana modes can be realized 
in one-dimensional systems has been inspired by the imminent possibility of discovering these fascinating objects. In 
conclusion, we believe that the studies presented here can be further developed in interesting theoretical directions as 
well as investigated experimentally. 



Acknowledgments 

For support, W.D. thanks the NSF under grant DMR 0644022-CAR and UChicago Argonnc, LLC, operator of 
Argonnc National Laboratory, under contract No. DE-AC02-06CH11357; M.T. thanks the CSIR, India; D.S. thanks 
DST, India under Project No. SR/S2/JCB-44/2010; S.V. thanks the Simons Foundation under Grant No.229047 and 
the U.S. Department of Energy, under Award No. DE-FG02-07ER46453. W.D. and S.V. thank the Indian Institute 
of Science for its generous hospitality at different stages of this work. We are grateful to Ilya Gruzberg, Shinsei Ryu 
and Jay Deep Sau for their illuminating comments. 



A. Kitaev, Pliysics- Uspekhi 44, 131 (2 001), |arXiv:cond-mat/0010440y 2 (2000). 

C. W. J. Beenakker. larXiv:1112.1950V 2 (20117 ' 

R. M. Lutchyn, J. D. Sau, S. Das Sarma, Phys. Rev. Lett. 105, 077001 (2010). 
Y. Oreg, G. Refael, and F. von Oppen, Phys. Rev. Lett. 105, 177002 (2010). 

L. Fidkowski, J. Alicea, N. H. Lindner, R. M. Lutchyn, and M. P. A. Fisher, Phys. Rev. B 85, 245121 (2012). 
A. C. Potter and P. A. Lee, Phys. Rev. Lett. 105, 227003 (2010). 

I. C. Fulga, F. Hassler, A. R. Akhmerov, and C. W. J. Beenakker, Phys. Rev. B 83, 155429 (2011). 

T. D. Stanescu, R. M. Lutchyn, and S. Das Sarma, Phys. Rev. B 84, 144522 (2011). 

S. Tewari and J. D. Sau, Phys. Rev. Lett. 109, 150408 (2012). 

M. Gibertini, F. Taddei, M. Polini, and R. Fazio, Phys. Rev. B 85, 144525 (2012). 

J. S. Lim, L. Serra, R. Lopez, and R. Aguado, Phys. Rev. B 86, 121103 (2012). 

M. Tezuka and N. Kawakami, Phys. Rev. B 85, 140508(R) (2012). 

R. Egger and K. Flensberg, Phys. Rev. B 85, 235462 (2012). 

S. Gangadharaiah, B. Braunecker, P. Simon, and D. Loss, Phys. Rev. Lett. 107, 036801 (2011). 

E. Sela, A. Altland, and A. Rosch, Phys. Rev. B 84, 085114 (2011). 

A. M. Lobos, R. M. Lutchyn, and S. Das Sarma, Phys. Rev. Lett. 109, 146403 (2012). 

R. M. Lutchyn and M. P. A. Fisher, Phys. Rev. B 84, 214528 (2011). 

A. M. Cook, M. M. Vazifeh, and M. Franz, Phys. Rev. B 86, 155431 (2012). 

F. L. Pedrocchi, S. Cliesi, S. Gangadharaiah, and D. Loss, Phys. Rev. B 86, 205412 (2012). 

D. Sticlet, C. Bena, and P. Simon, Phys. Rev. Lett. 108, 096802 (2012); D. ClievalUer, D. Sticlet, P. Simon, and G. Bena, 
Phys. Rev. B 85, 235307 (2012). 

P. San- Jose, E. Prada, and R. Aguado, Phys. Rev. Lett. 108, 257001 (2012); E. Prada, P. San- Jose, and R. Aguado, Phys. 
Rev. B 86, 180503 (2012). 

J. Khnovaja and D. Loss, Pliys. Rev. B 86, 085408 (2012). 

J. AUcea, Rep. Prog. Phys. 75 , 076501 (2012). 

T. D. Stanescu and S. Tewari. lMOCiv:1302. 54331 (2013). 

S. B. Ghung, H.-J. Zhang, X.-L. Qi, and S.-G. Zhang, Phys. Rev. B 84, 060510 (2011). 
V. Shivamoggi, G. Refael, and J. E. Mo ore. Phys. Rev. B 82, 041405(R) (2010). 
I. Adagideh, M. Wimmer, and A. Teker, arXiv:1302.26r2 (2013). 
J. D. Sau and S. Das Sarma, Nature Communications 3, 964 (2012). 

A. R. Akhmerov, J. P. Dahlhaus, F. Hassler, M. Wimmer, and C. W. J. Beenakker, Phys. Rev. Lett. 106, 057001 (2011). 
W. DeGottardi, D. Sen, and S. Vishveshwara, arXiv: 1208.0015 , to appear in Phys. Rev. Lett. (2013). 
W. DeGottardi, D. Sen, and S. Vishveshwara, New. J. Phys. 13, 065028 (2011). 
D. Sen and S. Vishveshwara, EPL 91, 66009 (2010). 

Y. Niu, S. B. Chung, C.-H. Hsu, I. Mandal, S. Raghu, and S. Chakravarty, Phys. Rev. B 85, 035110 (2012). 
J. D. Sau, C. H. Lin, H.-Y. Hui, and S. Das Sarma, Phys. Rev. Lett. 108, 067001 (2012). 
L.-J. Lang and S. Chen, Phys. Rev. B 86, 205135 (2012). 

P. W. Brouwer, M. Duckheim, A. Romito, and F. von Oppen, Phys. Rev. Lett. 107, 196804 (2011). 



30 



P. W. Brouwer, M. Duckheim, A. Romito, and F. von Oppen, Phys. Rev. B 84, 144526 (2011). 

V. Mourik, K. Zuo, S. M. Frolov, S. R. Plissard, E. P. A. M. Bakkers, and L. P. Kouwenhoven, Science 336, 1003 (2012). 
M. T. Deng, C. L. Yu, G. Y. Huang, M. Larsson, P. Caroff, and H. Q. Xu, Nano Lett. 12, 6414 (2012). 
L. P. Rokhinson, X. Liu, and J. K. Furdyna, Nature Phys. 8, 795 (2012). 

A. Das, Y. Ronen, Y. Most, Y. Oreg, M. Heiblum, and H. Shtrikman, Nature Phys. 8, 887 (2012). 

J. Alicea, Y. Oreg, G. Refael, F. von Oppen, M. P. A. Fisher, Nature Phys. 7, 412 (2011). 

S. Vishveshwara, Nature Phys. 7, 450 (2011). 

M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 (2010). 

X.-L. Qi and S.-C. Zhang, Rev. Mod. Phys. 83, 1057 (2011). 

C. Y. Teo and C. L. Kane, Phys. Rev. B 82, 115120 (2010). 
L. Fidkowski and A. Kitaev, Phys. Rev. B 83, 075103 (2011). 
K. T. Law and P. A. Lee, Phys. Rev. B 84, 081304(R) (2011). 
A. Altland and M. R. Zirnbauer, Phys. Rev. B 55, 1142 (1997). 

P. W. Brouwer, A. Furusaki, I. A. Gruzberg, and C. Mudry, Phys. Rev. Lett. 85, 1064 (2000); P. W. Brouwer, A. Furusaki, 
and C. Mudry, Phys. Rev. B 67, 014530 (2003). 

L A. Gruzberg, N. Read, and S. Vishveshwara, Phys. Rev. B 71, 245124 (2005). 

E. Lieb, T. Schultz, and D. Mattis, Ann. Phys. (NY) 16, 407 (1961). 

D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs, Phys. Rev. Lett. 49, 405 (1982). 
O. Motrunich, K. Damle, and D. A. Huse, Phys. Rev. B 63, 224204 (2001) . 

Q.-J. Tong, J.-H. An, J. Gong, H.-G. Luo, and C. H. Oh. larXiv:1211.2498l f2012). 
X.-G. Wen and A. Zee, Nucl. Phys. B 316, 641 (1989). 

A. Kitaev, Ann. Phys. (NY) 321, 2 (2006). 
D. S. Fisher, Phys. Rev. B 51, 6411 (1995). 

M. Ya. Azbel, Zh. Eksp. Teor. Fiz. 46, 929 (1964) [Sov. Phys. JETP 19, 634 (1964)]; D. R. Hofstadter, Phys. Rev. B 14, 
2239 (1976). 

Y. Last, Commun. Math. Phys 164, 421 (1993). 

S. Y. Jitomirskaya, Annals of Mathematics 150, 1159 (1999). 

F. Delyon, J. Phys. A 20, L21 (1987). 

G. Andre and S. Aubry, Ann. Israel Phys. Soc. 3, 133 (1980). 

B. Derrida and E. Gardner, J. Physique 45, 1283 (1984). 

F. M. Izrailev, A. A. Krokhin, and N. M. Makarov, Phys. Rep. 512, 125 (2012). 
P. J. Lloyd, J. Phys. C 2, 1717 (1969). 
D. J. Thouless, J. Phys. C 5, 77 (1972). 

C. Nayak, S. H. Simon, A. Stern, M. Freedman, S. Das Sarma. Rev. Mod. Phys. 80, 1083-1159 (2008). 



